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INTRODUCCION 


ESDE que el hombre diseñó y construyó el primer orde- 
nador siempre ha tenido el empeño de extender su uso a 
todas las áreas posibles de la vida humana. Pensemos en 
los primeros ordenadores, aun familiares para todos por 
las imágenes que de ellos hemos recibido a través de los 
medios de comunicación y espectáculos: nadie se queda- 
ba con la idea de la máquina de calcular que estaba en- 
cerrada en un edificio y que suministraba soluciones a las 
complicadas ecuaciones con las que los matemáticos la 
alimentaban a través de cientos de tarjetas perforadas. La imagen del or- 
denador era la del dispositivo que controlaba lo que le rodeaba, que sabía 
todo sobre todo y sobre todos nosotros, que era capaz de responder a cual- 
quier pregunta que se le formulara. Y esta imagen del ordenador alimentó 
nuestra mente y la de cientos de escritores de ciencia ficción que se dedi- 
caron a elucubrar sobre el tema. 

Literatura aparte, esta es la idea fundamental que debemos destacar al 
comenzar este libro. El ordenador es una máquina de propósito general, 
capaz de resolver problemas de la más diversa índole al usuario sin más 
que utilizar el programa conveniente y los circuitos de control adecuados. 
Hoy en día tenemos ejemplos claros en la mayoría de los aparatos que nos 
rodean: el ascensor automático que recuerda todos los pisos en los que 
debe parar, el horno que nos dice el tiempo de cocinado sin más que su- 
ministrarle el alimento a preparar y su peso, la máquina de coser con la 
que la abuelita hubiera soñado, que cose como ni a la más hábil costurera 
se le hubiera ocurrido hacerlo... 

En el campo de la ingeniería, el ordenador se supo aprovechar desde 
sus comienzos. Incluso cuando era demasiado voluminoso y demasiado 
caro para poder introducirlo en una lavadora, no lo era para las grandes 
empresas, para las universidades y centros de investigación. El técnico 
contaba con una nueva herramienta de trabajo. Todos los cálculos tedio- 
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sos, repetitivos que ocupaban la mayor parte de su tiempo podía descri- 
birlos en función de simples instrucciones. El ordenador las podía ejecu- 
tar tantas veces como él quisiera sobre distintos datos. Así se obtenían las 
soluciones que luego se podrían transformar en motores, alas de un avión, 
carenados o radioenlaces. 

Quizá movido a reflexionar por el mayor tiempo que el ordenador le 
permitía emplear en su labor creativa o por cualquier otra razón, el hecho 
es que cada vez el ingeniero le fue pidiendo más a su ordenador. Poco a 
poco lo fue programando para tareas más complejas. En especial, para las 
encaminadas a resolver una de las grandes preguntas que se plantea siem- 
pre el constructor de cualquier tipo de aparato: «¿Qué pasaría si...?» Así na- 
cieron los programas de simulación. 

¿Qué es simular? Según la Real Academia de la Lengua Española, si- 
mular es «Representar una cosa, fingiendo o imitando lo que no es». En- 
trando más en detalle en el tema que nos ocupa, podemos seguir la defi- 
nición que Shannon dio en 1975: «Simulación es el proceso de diseñar un 
modelo de un sistema real y realizar experimentos con este modelo con 
el propósito de comprender el funcionamiento del sistema o de evaluar di- 
ferentes estrategias para su operación.» 

Gracias a estas técnicas el ordenador nos puede mostrar lo arrugado 
que podría quedar nuestro flamante diseño de un avión-cohete si se ve ata- 
cado por una lluvia de meteoritos, o lo que sería de nuestro amplificador 
estereofónico si lo enchufáramos a una red de 380 voltios en vez de nues- 
tros habituales 220. 


La simulación tiene una importancia fundamental para el diseño en in- 
geniería. El diseño es el proceso por el cual, partiendo de unos requeri- 
mientos fundamentales o especificaciones, se obtiene un producto final que 
los satisface. Esto es general para todas las áreas de la ingeniería, y en par- 
ticular para la ingeniería electrónica. 


En el caso de la ingeniería electrónica, intentar abordar este diseño o 
síntesis de un circuito directamente es casi imposible en cuanto su com- 
plejidad sea mediana. Las ecuaciones a resolver nos llevan a unos sistemas 
complicadísimos en los que prácticamente nunca hallaremos una única so- 
lución para nuestras incógnitas (valores de componentes), aun suponien- 
do que pudiéramos llegar a ellas. Salvo en casos muy concretos, como son 
el diseño de filtros, no disponemos de métodos que permitan, a partir de 
unas especificaciones funcionales, llegar a la realización óptima (ni tan si- 
quiera aproximada) del circuito que las cumple. 


Es por esto que el diseño electrónico siempre se ha basado en métodos 
de ensayo/error/repetición. En la figura 1 se puede ver un sencillo orga- 
nigrama que muestra el proceso normal por el que se construye un cir- 
cuito de cualquier clase. Partiendo de unas especificaciones externas, da- 
das por el resto del equipo en el que debe ir integrado el circuito, por nor- 
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Fig. 1. Metodología de diseño de un circuito. 


mas internacionales o por las especificaciones comerciales del aparato a 
fabricar, se fija una arquitectura: 

Para esta arquitectura se hace un cálculo aproximado de los valores de 
los componentes a utilizar. Aquí, además de la teoría de circuitos, de las 
matemáticas, siempre ha jugado un papel importantísimo la experiencia 
del diseñador, llegando esta en algunos casos a sustituir totalmente a las 
otras. 

El siguiente paso es construir un prototipo del circuito, fácilmente mo- 
dificable. Sobre él se hacen las medidas necesarias para comprobar si cum- 
ple o no las especificaciones de partida. 

En caso de que no las cumpla, se debe estudiar la razón e introducir 
en el diseño las modificaciones necesarias para aproximarse a lo que se 
busca. 

Este proceso de medida/alteración se repetirá tantas veces como sea 
necesario hasta que nos encontremos dentro de los márgenes establecidos 
de funcionamiento, cuya rigidez dependerá del tipo de aplicación. El últi- 
mo prototipo servirá de base para la construcción de una preserie o del 
equipo final. 

Por este procedimiento se han desarrollado durante mucho tiempo los 
circuitos convencionales (placas de un televisor, receptores de radio...). 

Pero consideremos ahora otro tipo de diseños electrónicos: el de esos 
diminutos artefactos de color negro con muchas «patitas» por cada lado y 
que los expertos llaman «chips». Un chip o circuito integrado no es más 
que un circuito electrónico convencional miniaturizado hasta unos extre- 
mos que jamás hubiesen podido imaginar en los tiempos en que la válvula 
de vacío era el «best-seller» de los componentes. En un chip puede haber 
miles de transistores, ¡y más del 50% de la superficie que vemos es sólo en- 
voltorio! : 

Pensemos, pues, en el proceso de diseño de uno de estos circuitos. Sus 
características eléctricas le vienen dadas en gran parte por su tamaño. No 
se puede construir un prototipo con componentes convencionales y espe- 
rar que responda igual que el producto final. Construir un prototipo inte- 
grado es una labor muy costosa, pues lo que abarata tanto estos dispositi- 
vos es el gran volumen de fabricación. Incluso si lo pudiéramos construir, 
no sería posible realizar modificaciones sobre él, teniendo que realizar uno 
nuevo en caso de que existan errores en el diseño. 


La única salida posible en este caso es la simulación. Una simulación 
que considere todas las características especiales que inciden sobre este 
tipo de circuitos y nos ayude a respetar las muchas consideraciones de di- 
seño imprescindibles para realizar un dispositivo que funcione. 


Esta complejidad de los circuitos integrados ha sido el gran motor para 
el desarrollo de programas de simulación y la investigación sobre cómo 
utilizar el ordenador en el trabajo de concepción. 
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Para este fin es para lo que se han desarrollado los programas llama- 
dos de diseño asistido por ordenador / fabricación asistida por ordenador 
(más conocidos por sus siglas inglesas CAD / CAM : Computer Aided De- 
sign / Computer Aided Manufacturing). Prácticamente todas las áreas de 
la ingeniería en muy diversos campos en cada una de ellas poseen hoy pro- 
gramas de este tipo. 


En la figura 2 se puede ver de qué manera influye un ordenador en este 
proceso de diseño. La metodología es prácticamente la misma, pero tene- 
mos las siguientes ayudas: 


e La arquitectura a emplear muchas veces se reduce a elegir entre los 
distintos bloques funcionales de los que el ordenador ya posee en su me- 
moria de almacenamiento masivo. Por ejemplo, un bloque funcional pue- 
de ser «fuente de alimentación de 5 voltios y 10 amperios», que, en virtud 
de su uso en muy diversas situaciones, se puede haber construido con an- 
terioridad y emplear ahora directamente o con cambios muy pequeños. 

e Entrada de datos: Una vez elegida la arquitectura a emplear hay que 
construir los planos del circuito. El ordenador nos auxiliará en esta labor, 
permitiéndonos dibujar, de forma sencilla, los esquemas funcionales de 
los distintos elementos y ofreciéndonos procedimientos de edición del cir- 
cuito: conexionado, repetición de bloques, borrado de conexiones, inser- 
ción y eliminación de componentes... Todos estos datos los almacenará 
para utilizarlos posteriormente en las distintas partes del proceso. 


e Cálculo de componentes: No siempre es sencillo establecer cuál será 
su valor «a priori». El ordenador nos puede auxiliar con diversos métodos 
de cálculo para obtener una aproximación lo suficientemente buena para 
comenzar a trabajar. No obstante, es este un paso en el que la interven- 
ción directa del diseñador es imprescindible, pues el sentido común, en la 
mayoría de los casos, suele ser la mejor herramienta y la más cómoda de 
utilizar. 


+ Simulación: Esta es una de las principales ventajas que nos propor- 
ciona el diseño asistido por ordenador. En lugar de montar el circuito y 
realizar pruebas sobre él, proceso que puede resultar largo y costoso, el 
mismo ordenador nos dirá cómo responderá el circuito ante los estímulos 
externos que le fijemos. Para ello necesita saber cómo se comporta cada 
elemento básico. Esto lo conoce gracias a la biblioteca de modelos de dis- 
positivos de que dispone. El modelo que emplee, en especial para compo- 
nentes activos como transistores y diodos, es vital a la hora de obtener una 
simulación lo más próxima posible al caso real. 


e El proceso de corrección del circuito es ahora inmediato. Basta con 
acudir de nuevo al editor de entrada de componentes para alterar un va- 
lor de resistencia, cambiar de transistor o añadir un nuevo condensador 
de desacoplo. 
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Fig. 2. Proceso de diseño de un circuito asistido por ordenador. 


e Una vez que el análisis del circuito con el ordenador nos indica que 
las especificaciones se van a cumplir, debemos construir un prototipo para 
validar los cálculos. En caso de que el prototipo no responda a lo espera- 
do, debemos reconsiderar los modelos que se han empleado para la simu- 
lación y sustituirlos por otros más reales o más cercanos al caso. 


e Tras finalizar el diseño y validarlo con un prototipo, en el método tra- 
dicional quedaba una de las tareas más tediosas: el paso del diseño lógico 
a un diseño físico definitivo. La realización de una tarjeta de circuito im- 
preso es siempre una de las tareas más laboriosas en cualquier proceso. 


Hoy, partiendo del diseño lógico, que ya tenemos almacenado en nues- 
tro ordenador, podemos obtener el circuito impreso deseado de una for- 
ma casi automática. El ordenador nos auxiliará en las tareas de: 


e Colocación de los dispositivos de forma óptima. 


e Conexionado minimizando diversos parámetros: longitudes de pistas, 
tramos paralelos de señales, pasos de cara, problemas de ruido. 


e Generación de las máscaras para la realización física del circuito. 


En todo el proceso de evolución del trabajo de diseño podemos apre- 
ciar una característica fundamental, que muchas veces se nos olvida con- 
siderar a la hora de echar a volar nuestra imaginación: la mente del hom- 
bre está siempre por delante del trabajo que realiza el ordenador. No es 
el ordenador el que construye un yate, diseña un televisor o establece qué 
forma ha de tener el nuevo modelo de deportivo. El diseñador lleva a la 
práctica sus ideas utilizando distintas herramientas, y una de ellas y muy 
importante en la actualidad es el ordenador. 

En cualquier caso, conviene resaltar que normalmente nunca se habla 
de «diseño de circuitos por ordenador» simplemente, sino de «diseño de 
circuitos asistido por ordenador». En la palabra «asistido» se encierra gran 
parte del secreto de los ordenadores que, sólo en contadas ocasiones, es- 
capa del círculo de usuarios técnicos y se difunde en la sociedad. El dise- 
ño en ingeniería, no sólo electrónica, como ya hemos visto, normalmente 
tiene una cantidad de grados de libertad tan grande que el plantear un pro- 
blema general y hallar el método de programarlo en un ordenador es poco 
más o menos que imposible. Esto sólo se realiza en casos muy concretos 
y para pequeños problemas de los que se conoce muy bien la solución óp- 
tima y se puede delegar en el ordenador todas las tareas. En el resto de los 
casos, el diseñador se apoyará en el ordenador para ir realizando su tra- 
bajo. Este puede ayudarle en tareas de la más diversa índole: 


* Formas o módulos básicos que puede ensamblar, como un rompe- 
cabezas, para componer su producto. Para ello se ayudará de su gran ca- 
pacidad de almacenamiento para mantener una biblioteca de módulos. 
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e Normas de conexionado de módulos: puede detectar irregularidades 
en los conexionados que darán lugar a un mal funcionamiento. 

e Dibujo automático de los planos y guías de montaje, con la gran ven- 
taja de poder regenerarlos de nuevo tras un cambio sin producir el doble 
de trabajo de delineante. 

e Simulación del producto en las condiciones que se establezca, para 
evaluar si su comportamiento responde a las perspectivas para las que se 
diseñó. 

Así se podrían seguir dando ejemplos de situaciones en las que un or- 
denador puede asistir a la hora de diseñar cualquier cosa. 

El presente libro va a intentar introducirnos en las técnicas y herra- 
mientas destinadas a realizar las tareas de análisis de un circuito de modo 
automático, y a aprovechar los conocimientos que se obtienen del análisis 
para la ayuda en la síntesis. 

A lo largo de él se introducirán programas que intenten llevar a la prác- 
tica las ideas y los algoritmos matemáticos expuestos. El fin fundamental 
que persiguen es didáctico y, por tanto, en ellos se ha sacrificado en mu- 
chos casos la potencia o la rapidez de ejecución en favor de la claridad. 
Se pretende que el lector pueda extraer de ellos las ideas suficientes para 
atacar su propio problema de una forma más eficiente. 

Precisamente por esta razón es por lo que los programas se han de- 
sarrollado en lenguaje Pascal, fruto de los esfuerzos de la programación es- 
tructurada, y cuya validez como medio de transmisión de ideas en el cam- 
po de la programación está más que probada y avalada por su extendido 
uso, no sólo como lenguaje en sí, sino también como metalenguaje en des- 
cripción de algoritmos. Para facilitar la labor de comprensión al lector que 
no esté muy familiarizado con él se ha tratado de utilizar sólo las instruc- 
ciones más sencillas del repertorio disponible en cualquier buen compila- 
dor. De este modo se facilita también la tarea de adaptación para quien de- 
see construir sus futuros programas de diseño en otros lenguajes. 
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FUNDAMENTOS DEL DISEÑO DE 
CIRCUITO ASISTIDO POR 
ORDENADOR 


STE capítulo trata de uno de los temas que suele presen- 
tar menos interés al aficionado, pero cuyo conocimiento 
es imprescindible tanto para el profesional como para el 
que se quiere tomar un poco más en serio el pasatiempo 
de los ratos de ocio. Se trata, naturalmente, de las mate- 
máticas a emplear para poder representar un circuito de 
cara al ordenador y realizar cálculos sobre él. 

En todo este capítulo se supone que el lector posee 
unos conocimientos matemáticos sobre lo que es la deri- 
vación, los números complejos y las funciones trigonométricas, que no se 
desarrollan por quedar fuera del tema y la extensión del presente libro. 


E: 


DEFINICIONES FUNDAMENTALES 


Una red eléctrica es un conjunto de componentes tales como resisten- 
cias, condensadores, bobinas, etc., ligados entre sí mediante conductores. 
Una red es lineal cuando está constituida por elementos lineales y bi- 
laterales, que son en los que las ecuaciones que los representan son inde- 
pendientes de la dirección de la corriente o de lo grande que sea su valor. 
Elementos lineales típicos son las resistencias, bobinas y condensado- 
res (no electrolíticos), ya que en ellos no influye la posición en que se co- 
necten y se siguen comportando igual por muy grande que sea la corrien- 
te que los atraviese, siempre, por supuesto, ¡que no se lleguen a quemar! 
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Aunque más adelante se definirán formalmente los citados elementos, va- 
mos a ver aquí de una forma sencilla sus características fundamentales: 


e Una resistencia o resistor es un componente eléctrico que se «opo- 
ne» al paso de una corriente por él. Esta oposición es tanto mayor cuanto 
mayor sea su valor, que medimos con una unidad denominada ohmio. Una 
característica fundamental de estos componentes es que transforman par- 
te de la energía eléctrica que conducen en calor. 

e La bobina es un componente constituido por un conductor enrolla- 
do sobre un núcleo, normalmente cilíndrico. Su propiedad fundamental 
es la de almacenar energía en forma de campo magnético. 

e Un condensador es un componente constituido por dos conductores 
próximos entre sí y separados por un elemento aislante o dieléctrico. Al 
igual que la bobina, son capaces de almacenar energía, pero en este caso 
se realiza en forma de campo eléctrico. 


Estos circuitos los podemos ver representados en la figura 3. 

Elementos no linealés son aquéllos en los que no se dan las condicio- 
nes establecidas para los lineales. Ejemplos clásicos son los diodos y los 
transistores. 


e El diodo es un componente de dos terminales que permite que la 
corriente circule por él siempre que entre por el terminal que llamaremos 
ánodo y, por tanto, salga por su otro terminal, que llamaremos cátodo. En 
caso contrario, no deja pasar la corriente eléctrica. 

e Un transistor es un elemento de tres terminales. Hay distintos tipos 
de ellos, pero todos coinciden en que la corriente que circula entre dos 
de sus terminales se puede controlar con el tercer terminal: 


e En el caso de los llamados transistores bipolares, el control se realiza 
mediante la corriente que se hace pasar por el tercer terminal. Los 
terminales controlados se denominan colector (C) y emisor (E), y el 
terminal de control, base (B). Según los sentidos de las corrientes 
que lo atraviesan, hay dos tipos de transistores bipolares: 


e NPN: en ellos la corriente de control entra en la base, y la contro- 
lada circula de colector a emisor. 

e PNP: exactamente al contrario: la corriente sale de la base y la con- 
trolada circula de emisor a colector. 


e En los llamados transistores de efecto de campo (FET), es la tensión 
aplicada a uno de sus terminales, que llamamos puerta, la que con- 
trola la corriente entre los terminales llamados drenador (D) y fuen- 
te (S). 


Los circuitos a los que nos referiremos de momento serán lineales. Más 
adelante estudiaremos lo que se debe hacer para analizar los componen- 
tes no lineales. 
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Otros componentes básicos de una red son los generadores, elementos 
que introducen la señal en nuestro circuito. Los podemos dividir en dos 
grandes grupos: 


e Generadores independientes, cuyo valor de tensión no varía con el 
valor de cualquier otro parámetro de la red. 
e Generadores controlados, en caso contrario. 


En el siguiente apartado se estudian con más detalle dichos componen- 
es. 


_ 


v(t) E R v(t) L v(t) Cc 
v(t) =Ri(t) v(t)= L: —0— (y= Cc — A — 


Fig. 3. Elementos básicos de una red y ecuaciones que los representan. 


RELACIONES CORRIENTE-TENSION 


Tenemos como fin representar una red eléctrica de forma matemática. 
Debemos comenzar, por tanto, por estudiar la forma en que podemos re- 
presentar los elementos que la componen. Veamos la figura 3. En ella se 
ven los elementos más sencillos que podemos considerar. Por ellos circu- 
la una corriente que varía con el tiempo, que llamaremos i(t) y entre sus 
extremos hay una diferencia de tensión que denominaremos vít). Para de- 
finir su comportamiento eléctrico es necesario dar una ecuación que re- 
lacione estas dos magnitudes, denominada (¿se le ocurriría a usted un 
nombre mejor?) relación corriente-tensión. Así tenemos: 


e Resistencias: Su relación corriente-tensión viene dada por la Ley de 
Ohm en su forma más simple: la tensión en bornas de una resistencia es 
igual al producto de la corriente que por ella circula por una constante: 


v(t) = Ri(t) 


El valor de la constante es lo que llamamos resistencia, y se mide en 
ohmios. 
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e Bobinas: En este caso la relación no es tan simple. La tensión en bor- 
nas de una bobina es igual al producto de una constante por la derivada 
de la corriente que por ella pasa respecto al tiempo: 


di(t) 
dt 


La constante L es lo que se denomina valor inductivo de una bobina o au- 
toinducción, y su unidad de medida es el henrio. 

e Condensadores: En estos casos la relación es inversa a la de las bo- 
binas: la corriente que atraviesa un condensador es igual al producto de 
otra nueva constante por la derivada de la tensión en sus extremos respec- 
to del tiempo: 


v(t) = 


dv(t) 
dt 


En este caso la constante se denomina capacidad del condensador, y 
su unidad de medida es el faradio. 


i(t) =C 


+ 
v(t) i(t) 


Generador ideal de tensión Generador ideal de corriente 
Rg 
+ 
“(0 10) | Rg 
Generador real de tensión Generador real de corriente 


Fig. 4. Generadores de excitación de la red. 


Los elementos que actúan como excitadores de nuestra red reciben el 
nombre de generadores. En un amplificador el generador sería el disposi- 
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tivo que colocamos a su entrada: un micrófono o la salida del tocadiscos, 
por ejemplo. En la figura 4 se puede ver su representación sobre el papel. 

Un generador ideal es el que es capaz de suministrar su tensión o su 
corriente, según el tipo de que se trate, independientemente de las carac- 
terísticas y valores de los componentes que se coloquen a su salida. Estos 
generadores, como es lógico, no existen en la realidad. Sin embargo se es- 
tudian porque un generador real se puede considerar como un generador 
ideal con una determinada resistencia interna, que limita su funcionamien- 
to y lo hace dependiente del exterior. 


+ 
i(t) 
y eg(t) = K,:v(t) 
v(t) egít) eg(t) = K,:(t) 
+ 
(1) 
v(t) ig(t) = Kyv(t) 
ig(t) ig(t) = Kqi(t) 


Fig. 5. Generadores controlados ideales. 


Otros de los componentes normalmente empleados para representar 
circuitos electrónicos, sobre todo para el modelado de dispositivos com- 
plejos mediante elementos simples, son los generadores controlados. En 
la figura 5 se puede ver su diagrama y sus relaciones corriente-tensión. 
Como se puede apreciar, el valor de tensión o corriente que otorgan no 
es, como en el caso anterior, constante, sino que depende del valor de otra 
magnitud en cualquier parte del circuito. 
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También en este caso podríamos hablar de generadores controlados 
ideales y reales, al igual que para los generadores independientes. Sus re- 
presentaciones se ven en la siguiente figura: 


Rg 
+ 
10) 
+ 
MS eg(t) 


Fig. 6. Generadores controlados reales. 


Para terminar, veamos las ecuaciones que rigen el comportamiento de 
dos inductancias con acoplamiento mutuo o transformador. 


v2(t) 


L, L, 
Wi 
Fig. 7. Representación esquemática de un transformador. 
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Los puntos junto a las bobinas indican por cuál de sus extremos deben 
entrar las corrientes al transformador para que los flujos magnéticos crea- 
dos por ambos se sumen. 


Las relaciones corriente-tensión en este caso son de la forma siguiente: 
dil(t di2(t 
(0 y mM di2(0 
dt dt 


v2(t) = 12 di2(t) E: di1(t) 
dt dt 


v1(t) = L1 


Se pueden apreciar dos términos que relacionan ahora la corriente con 
la tensión en cada bobina: 


e El primero es el debido a su autoinducción, igual al que se tiene en 
caso de no existir acoplamiento. 

e El segundo es debido a la corriente por la bobina acoplada, y tiene 
la misma forma que el primero. El coeficiente M se denomina coeficiente 
de acoplamiento mutuo, y depende de la geometría del transformador. 


Este término representa la tensión que induce en un devanado la 
corriente que circula por el contrario. Esta tensión cambia de polaridad 
al cambiar el sentido de la corriente en el devanado inductor, lo que es 
importante a la hora de escribir las ecuaciones. Podemos establecer como 
norma general: 


Una corriente que entra por el extremo de la bobina marcado con 
el punto inducirá una tensión positiva en el extremo de la otra bobina 
también marcado con el punto, y viceversa. 


La siguiente figura nos muestran la diferencia de las ecuaciones en caso 
de que las corrientes por los devanados sean de distintos sentidos: 


+ + 
—+ AA 
i, i 
0) 
di di, 
be Ly MA 
Y; v 
di) di, 
va LM 77 
dt dt 
0) 
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7 

1 
di, diz 
a “a 
di> di, 
= "a 


vi(t) = L, 
v 


va(t) = L2 


+ + 
AS o 
1 

e ñ di di 

v; v> 
di di 

a 
1) 


Fig. 8. Flujos aditivos y sustractivos en un transformador. 


IMPEDANCIAS 


Como vemos, la ligadura entre las magnitudes que nos interesan no es 
simple. Aparecen derivadas. El análisis de un circuito por este camino nos 
llevaría a la resolución de un sistema de ecuaciones diferenciales, lo cual 
no es probablemente lo que más nos apetecía. 


v(t) , 


W= = o v(t) = V, : cos(wt) 


Fig. 9. Representación de una función sinusoidal en el tiempo. 
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Para evitar este conflicto demos un rodeo. Vamos a considerar que los 
generadores que atacan nuestro circuito siguen una variación de tipo si- 
nusoidal (ver figura 9). En este caso, y supuesto que ha pasado mucho tiem- 
po desde que se encendió el circuito (a esto se le llama régimen perma- 
nente), todas las corrientes y tensiones que por él circulen son también si- 
nusoidales. 


V, : cos(wt) 


Fig. 10. Representación de una función sinusoidal por un vector. 


Las señales sinusoidales se pueden representar por la proyección de un 
vector giratorio sobre un eje horizontal (ver figura 10). El vector lo pode- 
mos representar como un número complejo en forma módulo-argumento 
con la siguiente notación: 


vV ejot 


expresión en la cual ¡ es la unidad imaginaria (V -1). 
Esta forma de representar un número complejo es equivalente a otras 
que probablemente ya conocerá el lector, como pueden ser: 


e Forma binomial: a +j* b = cos(0wt) + j * sen(wt) 
e Forma módulo - argumento: V wt 


Conviene destacar también que es corriente en toda la teoría de circui- 
tos representar la unidad imaginaria con la letra j en lugar de con la letra 
i, como es habitual en matemáticas, debido a la confusión que originaría 
si se hiciera de esta última forma, por ser la letra i con la que normalmen- 
te se designan las intensidades de corriente. 

Según esta expresión, la señal sinusoidal representada en la figura 9 se- 
ría la parte real del vector así considerado: 


v(t) = Re [V ei*'] = Vo cos wt 


Normalmente en los cálculos no se toma esta parte real más que al fi- 
nal de las operaciones, y se trabaja durante todo el problema con las mag- 
nitudes complejas. 


ZE 


Si aplicamos ahora esta formulación a las relaciones corriente-tensión 
que se han visto para los componentes básicos, obtendríamos: 


V et =R]Ien para las resistencias 
5 jot 
V 60! = L A =j0L 1 elo para las bobinas 
¡ot 
I et = C OS =j0C V eje para los condensadores 


Si despejamos de estas ecuaciones el cociente entre los valores máxi- 
mos de tensión y corriente obtenemos: 


+ = = Zr en una resistencia 
V , ; 

ed =j0L =2ZIl en una bobina 

vV 1 

—=-— = Zc en un condensador 
I jot 


que son las expresiones de las impedancias de la resistencia, bobina y con- 
densador, respectivamente. Los cálculos se realizan para valores máximos 
de excitación, con lo cual no se maneja el término de la exponencial. En 
virtud de esto, las expresiones anteriores son las relaciones corriente ten- 
sión de los elementos básicos en régimen permanente sinusoidal. 

Es también interesante conocer que a los inversos de dichas cantida- 
des se les denomina admitancias. Así tendremos: 


1 ; ] ; : 
era = Yr Admitancia de la resistencia. 
r 
1 ; ; : 
Z = Y] Admitancia de la bobina. 
1 . ; 
ER = Yc Admitancia del condensador. 
Cc 
Z, Z, Z,=Z,+Z, 
AE A 
Serie 
Z; 
A 
2; Z, Z, 
Y, = Y, + Y, 
Paralelo 


Fig. 11. Asociación de impedancias. 
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Para las impedancias rigen las mismas normas de asociación que para 
las resistencias, esto es (ver figura 11): 


e Dos impedancias conectadas en serie suman su valor. 


e El inverso de la impedancia total de una conexión en paralelo de im- 
pedancias es igual a la suma de los inversos de los valores de las impedan- 
cias que se conectan. 


Fácilmente se ve, sin más que considerar su definición, que con las ad- 
mitancias estas reglas son a la inversa. 


Del mismo modo que transformamos las ecuaciones de los componen- 
tes sencillos, podemos transformar las de las inductancias acopladas, ob- 
teniendo unas ecuaciones de la forma: 


Vi =joL1 11 + joM 12 
V2 = joL2 12 + joM 11 


considerando los sentidos de corrientes y tensiones como en la figura 7. 


LEMAS DE KIRCHHOFF 


Ahora ya se han visto los componentes básicos de una red y las rela- 
ciones que los definen. Pasemos a estudiar ahora cómo atacar el proble- 
ma de su interconexión. 


Veamos un ejemplo sencillo de red: 


Rg 19 


ÉS 
E) 
p 
00000 
e 


Fig. 12. Ejemplo sencillo de una red lineal. 
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e Llamaremos malla a todo camino cerrado por el que pueda circular 
una corriente en nuestra red. En la red de la figura son mallas las formadas 
por los componentes: 


e Vg, Rg y L. 
e L,C yR2. 
e Vg, Rg, C y R2. 


Hay que llamar la atención sobre la última malla, normalmente no con- 
siderada, pero que tiene carácter de tal. 


e Llamaremos nudo al punto donde confluyen dos o más componen- 
tes. En la red de la figura son nudos los puntos de intersección de los com- 
ponentes: 


e Vg y Rg. 

e Rg,LyC. 
e CyR2. 

e Vg, L y R2. 


El último de ellos es un único nudo, aunque parezca «largo», debido a 
las líneas que unen los extremos aparentes de los componentes. 


Las leyes o lemas de Kirchhoff se refieren a estas dos entidades defini- 
das, y dicen así: 


1. La suma de las corrientes que entran en un nudo es igual a la suma 
de las corrientes que salen de él. 
2. La suma de las caídas de tensión en una malla es nula. 


Es muy importante señalar el sentido de corrientes y tensiones en cada 
componente a la hora de aplicar estas reglas. Normalmente se consideran 
las tensiones positivas en el mismo sentido que circula la corriente por el 
componente. 

Este sentido establecido «a priori» es arbitrario. Al resolver el proble- 
ma, si el sentido tomado era equivocado, el resultado obtenido será nega- 
tivo. 


Fig. 13. Asignación de corrientes a una red. 
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Para nuestro circuito, y considerando los sentidos de corrientes y ten- 
siones de la figura 13 tenemos: 


e Ecuaciones de corrientes: 


11 = 12 
12 = 13 + 14 
14 = I5 


e Ecuaciones de tensiones: 


Vg = Vrg + VI 
Vl = Ve + Vr2 


Se han omitido las ecuaciones del nudo Vg-L-R2 y de la malla Vg-Rg- 
C-R2 por ser redundantes con el resto, es decir, a partir de las presentadas 
se pueden obtener aquéllas sin más que realizar transformaciones algebrai- 
as. 

Con las ecuaciones establecidas y las relaciones corriente-tensión de 
cada componente se resuelve el circuito. En el caso más general, tomare- 
mos las relaciones corriente-tensión en su forma diferencial. La resolución 
del sistema de ecuaciones diferenciales al que llegaríamos nos daría la so- 
lución total del problema. Normalmente el problema no se ataca directa- 
mente así, sino que se divide en varios casos particulares que simplifican 
el estudio. 


o 


ANALISIS NODAL 


El método de resolución de circuitos más empleado en la mayoría de 
los programas de análisis por ordenador se deriva de la aplicación de los 
lemas de Kirchhoff de corrientes. Este es el método conocido como aná- 
lisis nodal, llamado así porque toma como variables del sistema de ecua- 
ciones a resolver las tensiones en todos los nudos que configuran el 
circuito. 

La única restricción necesaria para aplicar este método es que el cir- 
cuito no contenga fuentes ideales de tensión. En este caso se debe acudir 
a un análisis nodal modificado, que toma como variables independientes 
también las corrientes que dichos generadores ideales suministran al 
circuito. 

Esta restricción no es de vital importancia, ya que los generadores idea- 
les de tensión, por supuesto, no existen en la realidad, y en cualquier caso 
práctico vamos a tener generadores reales. 

Una vez establecida esta descripción, estableceremos aquí las normas 
prácticas para plantear las ecuaciones nodales de un circuito. El lector in- 
teresado en el proceso por el cual se obtienen estas reglas de manera teó- 


253 


rica lo puede encontrar en el «Análisis nodal de un circuito», en la 
página 121. 


En todo circuito debemos fijar en primer lugar un nudo de dato o de 
referencia, del cual no plantearemos ecuaciones, ya que si lo hiciéramos 
no introduciríamos información adicional para la resolución, pues podría- 
mos darnos cuenta de que dicha ecuación se puede obtener en función de 
las demás. Normalmente se toma el nudo marcado como «tierra» O «masa», 
por ser además la referencia natural del circuito, pero podríamos tomar 
cualquier otro. 


Para un circuito con n nudos (prescindiendo ya del de masa), el análi- 
sis nodal nos proporciona n ecuaciones con n incógnitas, que podemos re- 
presentar de la siguiente manera: 


Y11.V1 + Y12.V2 +... + Y1i. Vi +... + Yin. Vn = leql 
Y21.V1 + Y22.V2 +... + Y2i. Vi +... + Y2n. Vn = leg2 


AN 


Yn1 . V1 + Yn2 . V2 + .. + Yni . Vi +... + Ynn . Vn = legn 


o bien en forma reducida: 
[Yn] . [Vn] = [leg] 


Donde [Yn] es la matriz de admitancia nodal reducida, de dimensiones n 
por n; [Vn] es el vector de tensiones en los n nudos del circuito (nuestras 
incógnitas); e [leq] es el vector de generadores de corriente equivalentes 
que entregan corriente a cada uno de nuestros n nudos. 


e Determinación de los términos de la matriz de admitancias [Yn]. 


Cada uno de los elementos de esta matriz, que denominaremos gené- 
ricamente Yij, está constituido por una suma de términos, cada uno de los 
cuales representa la contribución al circuito de los elementos unidos a los 
nudos i, j o relacionados con ellos. 


En todos los casos se considera que las corrientes circulan partiendo 
del primer nudo citado en dirección al segundo, y que las tensiones son 
positivas en el primer nudo mencionado. Esto es: en un generador de 
corriente entre los nudos k y 1, la corriente irá del nudo k al l; si el gene- 
rador fuese de tensión, el nudo k sería el conectado al punto positivo de 
dicho generador. 
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1. Contribución de resistencias, bobinas y condensadores: Un elemen- 
to de este tipo, de admitancia Ya, conectado entre los nudos k y 1 del cir- 
cuito, da origen a cuatro términos en [Yn]: 


Fig. 14. Admitancias pasivas. 


Ykk: Ya Ykl: -Ya 
YIk: -Ya Yll: Ya 

2. Generador de corriente controlado por tensión: Supongamos que 
tenemos un generador de este tipo, de valor Gg, conectado entre los nu- 


dos k y 1, cuyo valor depende de la tensión existente entre los nudos r y s. 
Este generador introduce cuatro términos en [Yn], cuyos valores son: 


L 


K 
s 


Fig. 15. Generador de corriente controlado por tensión. 
Ykr: Gg Yks: -Gg 
Ylr: -Gg Yls: Gg 


3. Generador de corriente controlado por corriente: Un generador de 
este tipo, de valor Cg, conectado entre los nudos k y l, cuyo valor depende 


Ll 


de la corriente que recorre una admitancia de valor Yv conectada entre 
los nudos r y s, introduce cuatro términos en [Yn], de valores: 


L 


v Cg:i 


K 


Fig. 16. Generador de corriente controlado por corriente. 


Ykr: Cg. Yv Yks: -Cg . Yv 
Ylr: -Cg . Yv Yls: Cg.Yv 


4. Generador de tensión controlado por corriente: Dicho generador, 
de valor Dg y con una admitancia serie Yg, conectado entre los nudos k y 
l, cuyo valor depende de la corriente que recorre una admitancia de valor 
Yv conectada entre los nudos r y s, introduce cuatro términos en [Yn], de 
valores: 


Fig. 17. Generador de tensión controlado por corriente. 


Ykr: -Yg. Dv. Yv Yks: Yg. Dv. Yv 
Ylr: Yg. Dv. Yv Yls: -Yg. Dv. Yv 
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5. Generador de tensión controlado por tensión: En este caso, para 
un valor Bg con una admitancia serie Yg, conectado entre los nudos k y 
l, con un valor que depende de la corriente que recorre una admitancia 
de valor Yv conectada entre los nudos r y s, los términos introducidos en 
[Yn] tienen los siguientes valores: 


v 


Fig. 18. Generador de tensión controlado por tensión. 


Ykr: -Yg. Bv Yks: Yg. Bv 
Ylr: Yg. Bv Yls: -Yg. Bv 


e Determinación de los términos del vector de generadores de corrien- 
te equivalentes [leg]. 


Supondremos en este caso los mismos criterios de signos establecidos 
para la matriz de admitancias. 


1. Generador independiente de corriente, de valor Ig, conectado en- 
tre los nudos i y j: Da lugar a dos términos: 


J 


i 


Fig. 19. Generador independiente de corriente. 
legi: -Ig 
legj: Ig 
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2. Generador independiente de tensión, de valor Vg, conectado entre 
los nudos i y j, con admitancia en serie de valor Yg: Da lugar a dos 
términos: 


Yg 


Vg 


Fig. 20. Generador independiente de tensión. 
leqi: -Vg . Yg 
legj: Vg.Yg 


3. Generadores controlados: Introducen términos siempre que la 
rama que los controla contenga un generador de tensión independiente. 
Las siguientes figuras ilustran los distintos casos posibles: 


a) Generador de corriente controlado por tensión: 


Y: 
A j 
+ 
+ 
VE Vg ce Y, 
: i 
Fig. 21. Generador de corriente controlado por tensión. 
legi: -Gg . Vg 
legj: Gg.Vg 
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b) Generador de corriente controlado por corriente: 


Ey 
z j 
A 
Le 
- 
Vg C8virs 
i 
s 


Fig. 22. Generador de corriente controlado por corriente. 


+ 


leqi: -Cg . Vg . Yv 
legj: Cg.Vg.Yv 


c) Generador de tensión controlado por corriente: 


Ne 
r Vg 
j 
——» 
, 
+ + 
Vg Dg 1, 
i 
s 


Fig. 23. Generador de tensión controlado por corriente. 


legi: -Dg . Yg . Vg . Yv 
legj: Dg. Yg.Vg. Yv 
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d) Generador de tensión controlado por tensión: 


Y, Vg 


Fig. 24. Generador de tensión controlado por tensión. 


legi: -Bg . Yg . Vg 
legj: g.Yg. Vg 


PROGRAMA EJEMPLO 


A continuación se va a presentar un programa que realiza las siguien- 
tes funciones: 


1. Entrada de datos desde el teclado. 

2. Construcción de la matriz de admitancia nodal reducida y el vec- 
tor de generadores equivalentes. 

3. Presentación de ambos. 


Los tipos de componentes que acepta son: 


Resistencias 

Bobinas 

Condensadores 

Generadores independientes de corriente. 

Generadores independientes de tensión con admitancia interna resis- 


tiva. 

e Generadores de corriente controlados por la tensión o la corriente 
en una rama sin generadores. 

e Generadores de tensión con impedancia interna resistiva controla- 
dos por la corriente o la tensión en una rama sin generadores. 


Con estos componentes se pueden representar más del 90% de los cir- 
cuitos prácticos, por lo que no suponen una restricción muy grande. Se 
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pueden realizar fácilmente los cambios necesarios para soportar compo- 
nentes de otro tipo. 

El programa se halla profusamente comentado, con el fin de que el lec- 
tor pueda analizarlo sobre el mismo listado. No obstante, conviene aquí 
aclarar algunos puntos que pueden quedar oscuros: 


e La entrada de datos se realiza de forma interactiva, pidiéndose el pa- 
rámetro que en cada momento se debe dar para cada componente de for- 
ma explícita. 

e Los generadores de corriente se consideran con origen en el nodo 
del cual parte la corriente, y con destino en el nodo al que dicha corriente 
se suministra. 

e La matriz se presenta en forma compleja, con los términos inducti- 
vos separados de los capacitivos, para facilitar la labor a los programas de 
análisis que se presentarán en el capítulo siguiente, a la vez que se puede 
imprimir en una forma más compacta. 


El listado del programa es el siguiente: 


IS 


Número máxi 
UE 


(* El tercer indic: 


4 -a da iw —— [(Candensador) 
¡te Ccoericiente Cade Jw. _..—.. (LCOMndemsadot 


TAS naráme 
LOS parametr: 
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GGenerador : Real; 
(* Valor admitancia interna generador *) 
Error : Boolean; (* Mala entrada *) 


(* Procedimiento interno auxiliar para entrada de un 
generador de tensión *) 


procedure Gen Tension(var N1,N2: integer; 
var valor,G serie : real); 


begin 
writeln('Datos del generador:'); 
write('Terminal negativo conectado al nudo..>'); 
readin(N1); 
write('Terminal positivo conectado al nudo..>'); 
readln(N2); 


write('Valor del generador.....oooom........ == 
readln(valor); 

write('Valor del componente en serie........ 5 
readln(G serie); 

G_Serie := 1/G Serie; 


(* Por ser la entrada impedancia y requerir el 
análisis admitancias *) 
Valor:=Valor * G Serie; 
(* Generador equivalente de corriente *> 
end; 


(* Procedimiento interno auxiliar para entrada de un 
generador de corriente *) 


procedure Gen Corriente(var N1,N2: integer; 
var valor : real); 


begin 
writeln('Datos del generador:'); 
write('Terminal origen conectado al nudo....>'); 
readln(N1); 


write('Terminal destino conectado al nudo...>'); 
readln(N2); 
A O O AR a 
readlin(valor); 

end; 


(* Procedimiento interno y auxiliar para la entrada de 
un componente pasivo. Realiza también la entrada de 
los generadores controlados asociados a la rama que 
se introduce. *) 


procedure Entrada Componente (C: integer); 
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var 


Control : Char; 


begin 


TipoG:=0; 

Error:=False; 

writeln('Datos del componente:'):; 
write("Terminal positivo conectado al-nodo..>'); 
readin(N1)>; 

write(*Terminal negativo conectado al nodo..>'); 
readln(N2):; 

write('Valor del componente. oia. 3 
readin(ValorComponente); 


(* En caso de que el componente que se introduce sea 


una resistencia o una bobina se debe invertir su 
valor, pues la matriz que calculamos es de 
admitancias-*) 


if ((c= 0) or (c = -1)) then ValorComponente:= 
1/ValorComponente; 


(*=Se-procede ahora a la entrada de los parámetros del 
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generador controlado si es que existe *) 
Control:=*; 
while not-(Control- in PN*,*S*3)-do 
(*—Repetir-la pregunta mientras no se conteste "s" 
$ ue AY *) 
begin 


readin(Control); 
end; 
if Control ='S' then 
begin 


write('“Variable que lo controla (1/W= ===; 
readinítControl); 


case control of 
'v': TipoG:=1; 
ct 25 
else Error:=True; 
end; 


write("Tipo- de generador II eS 
readin(Control); 


case Control of 
Y: begin 


TipoG:=TipoG +-10; 
Gen—Tension(NC1,-NC2,-ValorGenerador; 
GGenerador); 
end; 
'1': begin 
TipoG:=TipoG-+-20; 
Gen—Corriente(NC1,-NC2,ValorGenerador); 
end; 
else 
Error:=True; 
end; 
end; 


(* Si no ha habido error en la entrada de datos se 
actualizan las matrices *) 


If£—not- Error then 
begin 
(* Actualización de la matriz de admitancias *) 
YINTNTCO= YN NEC VatorComponente; 
Y fN27N2C]:=Y[N2,N2,C]+ValorComponente; 
Y[NI,N2,C]J:=Y1N17N2,C3=ValorComponente; 
Y[N2,N1,C]:=Y[N2,N1,C]-ValorComponente; 


(*=Actualización de las matrices por la presencia de 
generadores controlados *) 


if TipoG-div-10-=-1-then 
begin 
(*=Generador de tensión. Lleva admitancia serie *) 
Y FNCIZNCIO]:+=YfNCI,NCI,O]+CGCGenerador; 
Y INC2,NC2,07:=Y[NC2,NC2,0]+CGenerador; 
Y[NC1,NC2,0]:=Y[NC1,NC2,0]-GGenerador; 
Y FENC27NC107:=YTNC2,NC1,0]-GGenerador; 
end; 
if TipoG mod 10 = 1 then 
begin 
(* Generador controlado por tensión *) 
Y[ÍNCI,N1,0]:=Y[NC1I,N1,O]J+ValorGenerador; 
Y NC2,N2,03+=YfNC2,N2,0]+ValorGenerador; 
Y[NC1,N2,0]:=Y[NC1,N2,0]-ValorGenerador; 
Y[NC2,N1,0]:=Y[NC2,N1,0]-ValorGenerador; 
end; 
if—TipoG-mod-10-=-2- then 
begin 
(* Generador controlado por corriente. Hay que 
multiplicar por-la admitancia del nudo 
controlador *;) 
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Y[NC1,N1,0]:=Y[NC1,N1,0]+ValorGenerador* 
ValorComponente; 
Y[NC2,N2,0]:=Y[NC2,N2,0]+ValorGenerador* 
ValorComponente; 
Y[NC1,N2,0]:=Y[NC1,N2,0]-ValorGenerador* 
ValorComponente; 
Y[NC2,N1,0]:=Y[NC2,N1,0]-ValorGenerador* 
ValorComponente; 
end; 

end 
else writeln('Entrada equivocada. ', 
'Repítala, por favor.'); 

end; 


(* Procedimiento de entrada de componentes *) 

begin 

(*- Se presenta un recordatorio de los tipos de 
componentes que se pueden introducir. Para acabar 


la entrada de componentes, pulsar la "gg" x) 


writeln('Tipos posibles de componentes: ' n= 


writeln(' - R: Resistencia. == 
writeln(' —- L: Bobina.'); 
writeln(' - C: Condensador.'); 
writeln('  - V: Generador indep. de tensión.'); 
writeln(' - I: Generador indep. de corriente.'); 
writeln(' - Q: Fin del circuito'); 
repeat 
WriteLn; 
Write('Introduzca el tipo de componente..... A 


ReadLn(Componente); 
(* Espera por el tipo de componente deseado *) 


(* Ahora se efectúan las operaciones necesarias para 
cada componente *) 


Case Componente of 


(* El tratamiento de R, L y C es análogo, y se realiza 
en el procedimiento "Entrada Componente"  *) 

'R' : Entrada _Componente(0); 

'L' : Entrada Componente(-1); 

!C* : Entrada Componente(1); 


(* Entrada de generador de tensión *) 
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(E Nisualizar la matriz de admitancias *) 


Procedure Ver Matriz Admitancias(var Y: Circuito); 
var-1,3]: integer; 


begin 
Writeln; 
Nriteln( "Matriz de admitancias del circuito”); 
Witt == 
Writeln; 
for i:=1 to NNodos do 
for-3:=1-to- NNodos- do- begin 


1 


ARA AA E RR E 
A RS 
end; 


end; 
(* Visualizar matriz de generadores independientes *) 
Procedure Ver Generadoresf[ var 1: Generadores); 


var 3j:integer; 


begin 
Writeln; 
Nriteln( "Matriz de generadores independientes”); 
MA A. == 
Writeln; 


for j:=1 to NNodos do 
A A AA 
end; 


(*=-FIN DE LAS RUTINAS AUXILIARES DE VISUALIZACION =*) 


( - __-_  -  -  _-  _  - _ _  _ _ => <Q€<»_E_IEOEOEE__E___EE_E_ A ) 


(*- Programa Principal *) 
BEGIN 


(*—Borrado de la matriz de admitancias *) 
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Para ilustrar el funcionamiento del programa se va a aplicar al cálculo 
de la matriz de admitancias reducida del siguiente circuito: 


Cc, 


Fig. 25. 
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En el circuito se han marcado los nodos con el número que se les va 
a dar al ejecutar el programa. Los valores asignados a los componentes son: 


Vag= 1V C1 = 10nF Ril= 1KQ 
Rg = 50Q C2 = 10pF R2 = 100Q 
gm = 0.3 L = 2mH 


A continuación se presenta el proceso tal y como se vería en la panta- 
lla del ordenador: 


Introduzca el número de nodos del 
circuito (Sin contar el de masa)...>4 


ATENCION: Ponga su teclado en mayúsculas 


Tipos posibles de componentes: 
- R: Resistencia. 


- L: Bobina. 
- C: Condensador. 
- V: Generador indep. de tensión. 
- I: Generador indep. de corriente. 
- Q: Fin del circuito 
Introduzca el tipo de componente..... >V 


Datos del generador: 
Terminal negativo conectado al nudo..>0 
Terminal positivo conectado al nudo..>1 


Valor del GENSrado Es: csi ale a >1 
Valor del componente en serie........ >50 
Introduzca el tipo de componente..... >C 


Datos del componente: 
Terminal positivo conectado al nodo..>1 
Terminal negativo conectado al nodo..>2 


valor, del ¡COMPONneNnte mi... ésos >10E-9 
Controla algún generador?............ >S 
Variable que lo controla (I/V)....... >V 
TÍpo de generador (I/V). 0... cr >1 
Datos del generador: 

Terminal origen conectado al nudo....>4 
Terminal destino conectado al nudo...>3 
Valor del ¡GENETAdO E coa. are ia a ar >0.3 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>2 
Terminal negativo conectado al nodo..>3 


Valor del :¡COMPOnenteS sii eicóimiaiaia aer es >1E3 
Controla algún generador?............ >N 
Introduzca el tipo de componente..... >L 
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Datos del componente: 


Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente.... 


Datos del componente: 


Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente.... 


Datos del componente: 


Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo.. 


v(LsL1) = 0.0000000000E+00 1/jw 


+ 2.0000000000E-02 + 1.0000001000E-01 


(TZ) = 0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + -1.0000000000E-08 


Y (1,3) = 0.0000000000E+00 1/j 


w 
E 0.0000000000E+00 + O.0000000000E+00 


Y(1,4) = 0.0000000000E+00 1/j 


wW 
+ 0.0000000000E+00 + -1.0000000000E-01 


Y(2,1) =  0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + -1.0000000000E-08 


Y(2,2) = 0.0000000000E+00 1/jw 


+ 1.0000000000E-03 + 1.0000000000E-08 


T(2,3) = 0.0000000000E+00 1/jw 


+  -1.0000000000E-03 + 0.0000000000E+00 


Y(2,4) = 0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + 0.0000000000E+00 


< 
a 
Y 
. 
dl 

Il 


0.0000000000E+00 1/jw 


+  -3.0000000000E-01 + 0.0000000000E+00 


Y(3,2) = 0.0000000000E+00 1/jw 


+ 2.9900000000E-01 + 0.0000000000E+00 


Y(3,3) = 5 .0000000000E+02 1/jw 
+ 1.0000000000E-03 + 0.0 


Y(3,2) = 0.0000000000E+00 1/jw 


+ 2.9900000000E-01 + 0.0000000000E+00 


Y(3,3) = 5 .0000000000E+02 1/jw 


+ 1.0000000000E-03 + 0.0000000000E+00 


jw 
jw 
jw 
jw 
jw 
jw 
jw 
jw 
jw 


jw 


jw 


jw 
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Y (3,4) 


+ O. 


Y(4,1) 
+ 3 
Y(4,2) 


+ -3. 


Y(4,3) 


me O: 


Y(4,4) 


+ 2. 


Lal 

_ 

19) 

— 
mun 


Como ya se especificaba en los comentarios del listado, cada elemento 
de la matriz de admitancias tiene tres componentes, que se deben al efec- 
to de las bobinas, resistencias y condensadores, respectivamente, que hay 
en el circuito. La matriz de generadores independientes es real. Esto im- 
plica que, en caso de haber varios en el circuito, todos tendrán la misma 
fase, que se tomará como origen de fases para el resto de corrientes y ten- 


siones. 
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.0000000000E-01 + 


= 0. 0000000000E+00 


0000000000E+00 + O. 


= 0.0000000000E+00 


= 0.0000000000E+00 


0000000000E-01 + O. 


0. 0000000000E+00 


0000000000E+00 + O. 


= 0.0000000000E+00 


0000000000E-02 + : E 


-1. 


1/jw 
0000000000E+00 
1/jw 
0000000000E-01 
1/jw 
0000000000E+00 
1/jw 
0000000000E+00 
1/jw 
0000000000E-01 


de generadores independientes 


2.0000000000E-02 
O. 0000000000E+00 
0.0000000000E+00 
0.0000000000E+00 


jw 
jw 
jw 
jw 


jw 


ANALISIS DE CIRCUITOS 
EN CONTINUA 


EGUN lo visto en el capítulo anterior, para analizar un cir- 
cuito es indispensable realizar el paso de los componen- 
tes físicos y su interconexionado a un sistema de ecuacio- 
nes. Esto se logra mediante 


e las relaciones corriente-tensión de los componen- 
tes lineales (o el modelo linealizado de los circuitos no li- 
neales), y 

e los lemas de Kirchhoff, que fijan las características 
de conexionado entre los primeros, 


resumidas mediante unas reglas para el análisis nodal, con el fin de im- 
plementarlas con facilidad por ordenador. 

Con esto se obtiene, en el caso más general, un sistema de ecuaciones 
diferenciales que representan el funcionamiento completo de la red. 

El programa de análisis descompone normalmente el problema en va- 
rios casos: 


e Análisis en continua. 
e Análisis en alterna. 
e Análisis transitorio. 


El presente capítulo va a abordar el primero de estos problemas, mien- 
tras que los dos restantes se tratarán en los siguientes. 
Se va a dividir en dos partes fundamentales: 


e Desarrollo de las matemáticas necesarias para llevar a cabo el análi- 
sis en corriente continua. 

e Programa sencillo para que el lector pueda comprobar una imple- 
mentación muy sencilla de los métodos matemáticos ilustrados en el capí- 
tulo. 
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EL PROBLEMA DE LOS COMPONENTES NO 
LINEALES 


El análisis de un circuito en continua es muy sencillo en el caso de cir- 
cuitos lineales. En este caso, las bobinas son un cortocircuito y los con- 
densadores un circuito abierto, con lo cual la red se simplifica enorme- 
mente. Solamente es necesario considerar resistencias y generadores. 

El mayor interés de este tipo de análisis está en el caso de circuitos no 
lineales. Es así porque, normalmente, los datos de funcionamiento que se 
obtengan para los componentes no lineales (punto de trabajo) van a fijar 
los parámetros para el circuito equivalente en alterna, donde normalmen- 
te nos interesan más los resultados. 

Un dispositivo no lineal es aquel cuya relación corriente-tensión sigue 
una variación no lineal, que se puede representar de la forma: 


v =f(i) 
Estos casos darían lugar a un sistema de ecuaciones diferenciales no li- 


neales difícil de resolver. Por esto lo que se hace es resolverlos por un pro- 
cedimiento iterativo de la siguiente forma: 


1. Se linealiza el circuito, aproximando la función por la tangente a 
ella en un determinado punto elegido arbitrariamente. Esto es: se toman 
como tensión y corriente unos valores que llamaremos 10, VO, tales que 
satisfagan que 


VO = f(10) 


Con ellos se aproxima la función por: 


v = VO + £(10) (i - 10) = Veq + Req i 


Veq = VO - f'(10).10 
Req = £'(10) 


que es la ecuación de un circuito como el de la siguiente figura: 
i 
—> Veqg 


+ v S 


Fig. 26. Circuito linealizado de un componente no lineal definido en tensión. 
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En caso de que el circuito no lineal nos viniese dado por una fórmula 
cuya variable independiente fuera la tensión y no resultase ser fácilmente 
despejable, el procedimiento sería análogo: 


i=f(v) 
10 = f (VO) 


¡ =10 + £(VO) (v - VO) = leq + — (1) 
Reg 


donde 


leq = 10 - £'(VO).VO 


Fig. 27. Circuito linealizado de un componente no lineal definido en corriente. 


2. Una vez linealizado se resuelve este circuito en continua. Con ello 
se obtienen unos valores de tensión y corriente para nuestro dispositivo, 
que no tienen por qué coincidir con los introducidos «a priori». 

3. Con los nuevos valores obtenidos se calcula el nuevo circuito equi- 
valente lineal del dispositivo. Este ciclo se repite hasta que la diferencia en- 
tre una solución y la siguiente sea menor que el margen de error que se 
pueda tolerar. 


A continuación se van a presentar las linealizaciones de algunos de los 
dispositivos no lineales más corrientes. 
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EQUIVALENTES LINEALES DE ALGUNOS 
DISPOSITIVOS 


Diodo de unión PN 


Sin entrar en efectos de segundo orden, el modelo más simple del dio- 
do que se puede establecer es el siguiente: 


+ 
a 


R id(t) = f(vd) = Is - [ev/Y - 1] 


SS 


Fig. 28. Diodo de unión PN y su relación corriente tensión. 


Las constantes de esta ecuación tienen el siguiente significado: 


e Is: Corriente de saturación. La podemos hallar en las características 
del catálogo, por representar la corriente que circula por el diodo bajo po- 
larización inversa. 

e Vt = KT/q , en donde: 


K: Constante de Boltzmann: 1.38? x 10 Jul/Kelvin 
T: Temperatura en Kelvin 
q: Carga del electrón: 1.6 x 10"? Q 


Con todo esto, a temperatura ambiente, 
Vt = 0.025 V. 

Junto con la ecuación del diodo tenemos la ecuación obtenida por la 

ley de Kirchhoff de tensiones: 
Vg = vd + id. R (3) 

El procedimiento más corriente consiste en realizar una construcción 
gráfica (ver figura 29). Sobre la curva que representa el comportamiento 
del diodo se traza la recta que fijan las condiciones externas al diodo (cuya 


ecuación nos ha dado la ley de Kirchhoff), llamada recta de carga. El pun- 
to de intersección de ambas es la solución del problema. 
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i¡=1, [e*/* -1] 


Ids 


Vds Vg 


Fig. 29. Resolución gráfica. 


Para emplear el procedimiento iterativo expuesto en este capítulo se 
debe linealizar la ecuación del componente. Su derivada es: 


Fw) = == (evM- 1) 


Tomemos ahora una solución arbitraria que denominaremos IdO0, VdO, 
tales que verificarán que: 


IdO = Is (ev40/M: — 1) 


Particularizando la derivada en este punto se obtiene: 
f'(VO) = e 3 (ev0M - 1) = gm 
vt 


La ecuación linealizada del diodo será, según (1): 
id = Ido + gm ( vd VdO ) (4) 
Sustituyendo el valor de id que se obtiene de esta ecuación en la (3), 
se resuelve el problema en esta primera aproximación: 


Vg = (Id0 + gm (vd Vd0)) R + vd 
= (Id0 gm VdO) R + (gm + 1) vd 


Vg — (IdO — gm Vd0) R 
gmR+1 


vd = = Vd1 
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Y sustituyendo este valor de vd en (2): 
Id1 = Is (evi - 1) 


se obtiene una nueva solución aproximada Id1, Vd1, con la cual se repe- 
tiría el proceso. 


Transistor bipolar 


El modelo más sencillo para el transistor es el de Ebers-Moll, que nos 
proporciona un circuito equivalente como el de la figura: 


Fig. 30. Circuito equivalente del transistor bipolar NPN. 


Este es el circuito general, empleado para el análisis transitorio, y en 
ese capítulo se volverá sobre él. Aquí, de momento, los condensadores se 
eliminarán por ser un circuito abierto en continua. 

Aplicando la ley de Kirchhoff de corrientes en los tres nudos Emisor 
(E), Colector (C) y Base (B), junto con la ecuación del diodo, citada ante- 


riormente, se obtiene: 
ie = ide + Ar . idc = 
= les (eVbeM — 1) + ar los (eveeM — 1) 
ic = — idc + af. ide = 
= — Ics (evbeM — 1) + af les (evbe/M: — 1) 
ib = ic —- ie 
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Valores típicos para estos parámetros son: 


les = Ics = 1 nA 
af =.98 
ar =.5 


Entre af y la ganancia en corriente del transistor, f, existe la relación: 


EA 
B+1 


En este caso, al haber dos ramas con elementos no lineales, la resolu- 
ción se complica bastante. 

Tomaremos para ilustrar el método a seguir un sencillo circuito con 
un transistor en configuración de emisor común: 


af 


Fig. 31. Circuito bajo análisis y circuito equivalente. 


DL 


Aplicando la ley de Kirchhoff de corrientes a los nudos de base y co- 
lector, se obtiene: 


Ib + af ide + ar idc — ide - idc = 0 
lc + idc — af ide = 0 


Operando con esas ecuaciones resulta: 


(Vgb - Vb) (1/Rb) + (af - 1) ide + (ar- 1) idc=0 (5) 
(Vgc — Ve) (1/Rc) af ide + idc = 0 (6) 


Se toma una solución arbitraria para las tensiones base-colector y ba- 
se-emisor Vbe0 y Vbc0. Aproximando por la tangente a la curva caracte- 
rística de los diodos base-emisor y base-colector, como se hizo en (4), obte- 
nemos: 


ide = Ide0 + gme (vbe - Vbe0) (7) 
gme = (les/Vt) (eve/M — 1) 
idc = IdcO + gmc (vbc — VbcO) (8) 


gmc = (Ics/Vt) (eveM — 1) 


El circuito equivalente lineal para esta aproximación será: 


Fig. 32. Circuito linealizado en la primera aproximación. 
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En el cual se ha hecho: 
lege = Ide0 -— gme Vbe0 
legc = Idc0 - gme Vbco 


Este circuito lineal se puede resolver por el método del análisis nodal, 
y obtener unos valores de Vb y Vc. Con ellos podremos calcular la segun- 
da aproximación para las tensiones vbe y vbc, ya que, por construcción 
del circuito: 
vbe = Vb 
vbc = Vb — Ve 


Con esta segunda aproximación se repetiría de nuevo el cálculo completo. 


Transistor de efecto de campo 


Su modelo más corriente se puede representar: 


G 


Fig. 33. Modelo en continua del transistor FET. 
El valor de la corriente que circula entre drenador y fuente, que lla- 
maremos Id es el siguiente: 


Id Idss == ( 51 A. EA 1) Para la zona lineal. 
Pp 


Id = Idss (1 — E y Para saturación. 
Pp 


En esta expresión: 


e Vp es la tensión de estrangulamiento o de «Pinch-off». Es el valor de 
la tensión puerta-fuente por encima del cual la corriente drenador-fuente 
no aumenta. 


e Idss es la corriente de drenador para Vgs = 0. 


IS 


Para el análisis en continua las capacidades serán un circuito abierto, 
con lo que desaparecerán del modelo, simplificando así el cálculo. 


A 


Transistor MOSFET 


Se expone aquí el modelo de primer orden, o de Shichman Hodges: 


Fig. 34. Modelo en continua del transistor MOSFET. 


La fórmula para ID es la misma que en el transistor FET de unión. 

Del mismo modo que en el transistor FET, las capacidades en continua 
desaparecerán del modelo, simplificando así el cálculo. Vemos como, al 
desaparecer las capacidades, nos queda la puerta (G) sin conectar a nin- 
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gún nodo. Por esto es por lo que a los transistores MOSFET se les deno- 
mina transistores de efecto de campo de puerta aislada. 


RESUMEN 


El método aquí expuesto puede no llevarnos a la solución en todos los 
casos que se nos ocurra plantear. Como todo método iterativo, su conver- 
gencia hacia la solución del problema sólo está garantizada bajo determi- 
nadas condiciones, en las que aquí no entraremos. Baste con decir que, si 
la solución inicial es lo bastante próxima a la real, el método suele dar re- 
sultados muy buenos y en pocas iteraciones. Pero, en caso contrario, pue- 
de incluso diverger, alejándonos cada vez más de nuestro objetivo. 

En estos casos se emplean otros métodos, como son el método de New- 
ton modificado y el método de Broyden, para cuya consulta se remite al 
lector interesado a la bibliografía básica que podrá encontrar al final del 
libro. 

Otro problema que se puede plantear es el de Overflow, esto es, encon- 
trar un valor de tensión o de corriente que desborde la capacidad de nues- 
tro ordenador. Esto puede ocurrir, por ejemplo, si se llega a una iteración 
que presenta un valor demasiado alto para la tensión en un diodo o en la 
base de un transistor. Este problema se puede resolver fácilmente si se to- 
man unos límites superiores para estos valores, fijados a partir de las con- 
diciones que imponga el circuito en el peor caso. 


PROGRAMA EJEMPLO 


Para ilustrar los métodos expuestos en este capítulo se presenta un pro- 
grama capaz de analizar circuitos lineales en continua y circuitos con dio- 
dos. 

Este programa incorpora para el análisis los siguientes procedimientos: 


e Entrada de un circuito y visualización de las matrices que lo definen. 
Este procedimiento es el mismo que se presentó en el capítulo anterior, 
en el que se ha suprimido la parte de entrada de inductancias y conden- 
sadores, no necesaria en este caso, con lo cual la matriz de admitancias in- 
terna se reduce en una dimensión. 

e Resolución de un sistema de ecuaciones por el método de descom- 
posición L-U. Estos procedimientos se detallan en el capítulo «Resolución 
de sistemas de ecuaciones», donde forman parte integrante de los progra- 
mas que allí se estudian. 


Ss) 


e Procedimientos para la entrada de los diodos, linealización y análisis 
del circuito. Son los propios del capítulo presente, y, por ello, se presen- 
tan con un gran número de comentarios para su fácil seguimiento. Por el 
contrario, en los programas ya tratados éstos se han suprimido, para con- 
centrar la atención en la parte original que aquí se pretende analizar. 


Conviene aclarar, al margen de los comentarios del programa, los si- 
guientes aspectos: 


e Para la tensión en los diodos se han tomado unos límites máximo y 
mínimo. Debido a la característica exponencial del diodo, para tensiones 
grandes, el valor que toma la función exponencial involucrada en el 
cálculo es enorme, pudiendo producir «Overflow» en nuestro ordenador. 
Del mismo modo puede existir «Underflow» para tensiones negativas, pues 
el valor de la exponencial es muy pequeño. 


Por estas razones se han tomado unos límites superior e inferior con 
los cuales se truncan dichos cálculos. Estos valores dependerán de la ca- 
pacidad de manejo de datos del compilador que estemos usando. 


e Las rutinas de entrada se han mantenido con una estructura básica 
igual a la expuesta en el capítulo «Fundamentos del diseño de circuitos 
asistido por ordenador», en la página 13. En este caso se usa sólo la parte 
correspondiente a los términos reales, desapareciendo, por tanto, la ter- 
cera dimensión de la matriz que indicaba el tipo de componente introdu- 
cido. El usuario debe tener esto en cuenta, sustituyendo en el circuito que 
desee analizar las bobinas por cortocircuitos y los condensadores por cir- 
cuitos abiertos, realizando así la asignación de nudos. 


e Los valores que se asignen a los diodos como condiciones iniciales 
son sumamente críticos. No se deben asignar tensiones nulas, por produ- 
cir errores en la resolución del sistema de ecuaciones. Valores en torno a 
0.4 voltios pueden considerarse recomendables, si bien pueden necesitar- 
se varios intentos para llegar a obtener un resultado aceptable. 


Una vez establecidas estas consideraciones, he aquí el programa: 


== rasrriceicnes 
A 


(* Constantes *) 
Const 


Max_Nodos =- 20; 
(+*=Número máximo de nodos del circuito *) 


(* Tipos del sistema *) 


Type 


(* Tipo de matriz que contendrá las admitancias del 
circuito. En este caso es del mismo tipo que-la 
matriz del sistema de ecuaciones *) 


Sistema = Array [0..Max_Nodos,0..Max-_Nodos] 
of real; 


(* Tipo de matriz para los generadores equivalentes *) 


Vector = Array[0..Max_ Nodos] of real; 


(* Procedimiento para la entrada del circuito 
Los parámetros que usa son: 
X: Matriz de admitancias del circuito: 
Parámetro de salida del procedimiento- 
Y: Matriz de generadores del circuito. 
Parámetro de salida del procedimiento. 
N: Número de nodos del circuito: 
Parámetro de entrada. 


=: 


(* Procedimiento para la entrada del circuito 
Los- parámetros que usa son: 
X: Matriz de admitancias del circuito: 
Parámetro de salida del procedimiento. 
Y: Matriz de generadores del circuito. 
Parámetro de salida del procedimiento. 
N: Número de nodos del circuito. 
Parámetro de entrada. 
3 


Procedure Entrada_del Circuito ( var Y: Sistema; 
var I: Vector); 
var 
Componente ==> (*—Tipo de componente *,) 


SZ 


Ni, (* Nodo origen del 


N2, (* Nodo destino del 
NC1, (* Nodo origen del elemento 
NC2, (* Nodo destino del elemento 
TipoG :—Integer; 

(* Tipo Generador 
ValorComponente, (* Valor componente 
ValorGenerador, (* Valor generador 
GGenerador : Real; 


componente *) 
componente *) 
controlado *) 
controlado *) 


controlado *) 
en la rama *) 
controlado *) 


(* Valor admitancia interna generador *) 
Error : Boolean; (* Mala entrada *) 


(* Procedimiento interno auxiliar para entrada de un 


procedure Gen Tension(var N1,N2: integer; 


generador de tensión *) 


var valor,G serie 


begin 


writeln('Datos del generador:'); 


real); 


write('Terminal negativo conectado al nudo..>'); 


readlin(Nl1); 


write('Terminal positivo conectado al nudo..>'); 


readln(N2); 


write('Valor del generador............. 


readin(valor); 


write('Valor del componente en serie... 


readln(G_serie); 
G_Serie := 1/G Serie; 


(* Por ser la entrada impedancia y requerir el 
análisis admitancias *) 


Valor:=Valor * G Serie; 


(* Generador equivalente de corriente *) 


end; 


(* Procedimiento interno auxiliar para entrada de un 
generador de corriente *) 


procedure Gen Corriente(var N1,N2: integer; 
var valor : real); 


begin 
writeln('Datos del generador:'); 
write('Terminal origen conectado al nudo....>'); 
readin(Nl1); 
write('Terminal destino conectado al nudo...>'); 
readln(N2); 
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write Yator-del generador == 23 
readlin(valor); 
end; 


(* Procedimiento interno y auxiliar para la entrada de 
una resistencia. Realiza también la entrada de los 
generadores controlados asociados a la rama que se 
introduce.  *) 


procedure Entrada Componente; 


var 
Control char; 


begin 
TipoG:=0; 
Error:=False; 
writeln('Datos del componente:'); 
write('Terminal positivo conectado al nodo..>'); 
readlin(Ni1); 
write('Terminal negativo conectado al nodo..>'); 
readln(N2); 
write('Valor del componente................. >); 
readin(ValorComponente); 


(* Se invierte el valor de la resistencia por estar 
calculando una matriz de admitancias *) 


ValorComponente:=1/ValorComponente; 


(* Se procede ahora a la entrada de los parámetros del 
generador controlado si es que existe *) 
Control:=' '; 
while not (Control in ['N','S']) do 
(* Repetir la pregunta mientras no se conteste "s" 
o "N" *) 
begin 
write('Controla algún generador?............ === 
readln(Control):; 
end; 
if Control ='S' then 
begin 
write('Variable que lo controla (1/V)....... => 
readin(Control); 


case control of 
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Nx PipoG==1; 

=p 727 

else Error:=True; 
end; 


write("Tipo-de-generador-(1/1) 5 =$ 
readln(Control); 
case Control of 
tv; begin 
TipoG:=TipoG + 10; 
Gen—Tension(fNCi,—NC2, ValorGenerador, 
GGenerador>); 
end; 
'*1'; begin 
TipoG:=TipoG + 20; 
Gen Corriente(NC1,-NC2,ValorGenerador); 
end; 
else 
Error:=True; 
end; 
end; 


(* si no ha habido error en la entrada de datos se 
actualizan las matrices *) 


If—not-Error-—then 
begin 
(* Actualización de la matriz de admitancias *) 
YN N1-=Y[NI,NiHValorComponente; 
Y[N27,N2]:=Y0N2,N2]+ValtorComponente; 
YFNI7N237=Y[N17N23=ValorComponente; 
Y[N2,N1]:=Y[N2,N1]-ValorComponente; 


f*Actualización de las matrices por la presencia de 
generadores controlados-*) 


if TipoG div 10 =- 1- then 
begin 
(*—Cenerador-de tensión. Lleva —admitancia—serie-*) 
Y FNCI7NCI=YfFNCI/NCIP*CGenerador; 
Y ÍNC2,NC2]:=Y[NC2,NC2]+GGenerador; 
Y[NC1,NC2]:=Y[NC1,NC2]-GGenerador; 
YTNC2,NC1]==YINC2,NC1]-GGenerador; 
end; 
if TipoG mod 10 = 1 then 
begin 
(*_Generador controlado por tensión *) 
Y FNCI7;N1J:=Y ENCINA ValorGenerador; 
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Y[NC2,N2]:=Y[NC2,N2J+ValorGenerador; 
Y ÍNCI7N2]+=Y[FNC1,N2]=ValorGenerador; 
Y[NC2,N1]:=Y[NC2,N1]-ValorGenerador; 
end; 
if TipoG-mod-30-=-2-—then 
begin 
(* Generador controlado por corriente. Hay que 
multiplicar por la admitancia del nudo 
controlador *) 
Y[NCI,N1]:=Y[FNCI;Ni]+ValorGenerador* 
ValorComponente; 
Y[NC2,N2]:=Y[NC2,N2]+ValorGenerador* 
ValorComponente; 
Y[NC1,N2]:=Y[NC1,N2]-ValorGenerador* 
ValorComponente; 
Y FNC2,NIJ3=YfNC2,N1]=Valorcenerador* 
ValorComponente; 
end; 
end 
else witeln( Entrada equivocada.” 
"Repítala, porfavor. ); 
end; 


(* Comienzo del procedimiento entrada de componentes*) 


begin 


(*- Se presenta un recordatorio de los tipos de 
componentes que se pueden introducir. Para acabar 
la entrada de componentes, pulsar la "o" x) 


writeln('Tipos posibles de componentes:'); 
writeln(' “- R: Resistencia.'):; 
writelnf—=1:Generador-indep- de tensión); 
writeln(-=-1:5-Generador-indep. de corriente.); 
writeln(* ——- Q: Fin del circuito'); 


repeat 
WriteLn; 
Write( "Introduzca el tipo de componente... >) 
ReadLn(Componente); 
(Espera por el tipo de-componente- deseado-—*) 


(*Ahora-se efectúan las operaciones necesarias para 
cada componente *) 


Case Componente of 
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(* El tratamiento de R, L y C es análogo, y se realiza 
en el procedimiento "Entrada Componente"  *) 
"R'-: Entrada Componente; 


(* Entrada de generador de tensión *) 
'v' : begin 

Gen Tension(NC1, NC2, ValorGenerador, 
GGenerador); 

(* Introducción de la admitancia serie *) 

Y[NC1,NC1]:=Y[NC1,NC1]+GGenerador; 

Y[NC2,NC2]:=Y[NC2,NC2]+GGenerador; 

Y[NC1,NC2]:=Y[NC1,NC2]-GGenerador; 

Y[NC2,NC1]:=Y[NC2,NC1]=-CGenerador; 

(* Efecto en el vector de generadores *) 

I[NC1]:=1[NC1]-ValorGenerador; 

1[NC2]:=1[NC2]+ValorGenerador; 

end; 


(* Entrada de generador de corriente *) 
115 begin 


Gen Corriente(NC1, NC2, ValorGenerador); 
I[NC1]:=I[NC1]-ValorGenerador; 
I[NC2]:=I[NC2)]+ValorGenerador; 
end; 
end; 
until Componente='0Q*' 
end; 


Procedure Resuelve Sistema(var A: Sistema; 
var B: vector; Orden: Integer;var X:Vector); 


var 
i,j,k : integer; 
E : Vector; 


(* Solución intermedia de LY = B *) 
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A A e A A e 
cientes matriz Il, ”* 


Irsa. >olucion Y 


Const 


Ids =-1E=-9; 
(*—Corriente de saturación de los diodos 
Max_Diodos = 20; 
(* Número máximo de diodos en la red 


VMax =-2;-(* Valor máximo tensión en diodos * 
Vmin ==13(* Valor mínimo tensión en diodos 
Type 


(*— Descripción diodos en el circuito 
Diodos = Array([1l..Max_Diodos] of 


record 
Na, (* Nodo del ánodo *) 
Nc:—iínteger; (*_Nodo del cátodo *) 
Va: Real; (* Tensión en la iteración actual *) 
end; 


Function id(V:- real): real; 
(+ Evaluación de la corriente por el diodo 
Begin 
If V>Vmax then V:=Vmax; 
(* Límite máximo para evitar el overflow 
If V<Vmin 
then Id :=-= Ids 
(4 Límite mínimo para evitar el under£low 
else Id:=-Ids * (exp (-V-/-0,025-3)--1); 
End; 


Function DId(VO: real): real; 
(* Evaluación de la derivada del diodo en V=VO 
Begin 
TI£—VO>Vmax then VO:=VWmax; 
If VO<Vmin 


then VO := Vmin; 
(* Límite mínimo para evitar el underflow 
Did:=-Tds-—*——exp-(-V0-/-0.025-3)-/0.025;5 
End; 


(Procedimiento para la entrada de los diodos *) 


Procedure Entrada Diodos(var D: Diodos; 
var ND: Integer); 
Var 
i-: integer; 
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id: 


== 


==: 


=. 


Begin 


For-1:=1 to ND-do 
begin 

WriteLn; 
WriteLn(*Diodo 1,1: 5); 
E O A == 
ReadEn(D[fi].Na); 
A A A === 
ReadLn(D(i].Nc); 
repeat 


Write('Tensión primera aproximación. ........ >); 


Readin(DHiTVa); 
until Abs(D[(i].Va) > 1E-10; 
end; 
end; 


(*—Procedimiento para presentar los resultados *) 


Procedure Presenta_Solucion (var V: Vector; 
N: Integer); 

var 

i:—integer; 
begin 

For i:= 1 to N do 

Writeln('V',i,' = ',V([i]); 

end; 


(* Procedimiento que realiza el análisis del circuito 
no lineal. Parámetros: 

E : Matriz de admitancias reducida: 

1 Vector de generadores independientes. 

NN : Número de nodos. 

D : Matriz de descripción de los diodos. 

ND : Número de Diodos. 

E Error-máximo permisible. 

= MI += Número máximo de iteraciones. *) 


Procedure Analisisftvar Y: Sistema; var G:Vector; 
var-NN: Integer; var D: Diodos; 
ND: Integer; E: real; 


MI: Integer); 


Var 

(4 Matriz real con los coeficientes del sistema de 
ecuaciones *) 
A : Sistema; 


B, (* Términos independientes *) 


Vi : Vector; (* Solución iteración i *) 
NIt : Integer; “(* Número iteración actual *) 
Final : Boolean; (* Indica fin de la iteración *) 
Iter  : Boolean; 
i,j,k : integer; 

F,DF —: Real; (* Variables auxiliares *) 

Begin 

Final := False; 

Nit =1; 

If ND= 0 


(*-Si no hay diodos el circuito es lineal *) 
then begin 
(*_Se copia el sistema de partida *) 


Ar=Y; 
B:=G; (* Términos independientes *) 


(* Resolución del sistema en la presente iteración *) 
Resuelve sistema(A,B,NN,Vi); 
(* Se presentan los resultados *) 
Writeln; 
Writeln('Solución del circuito:”'); 
Presenta Solucion (Vi,NN); 


end (Fin de tratamiento del análsis lineal *) 


else begin 
(*—Hay diodos. El análisis es no lineal *) 


while ((NIt<=MI) and (not Final)) do 
begin 


(*-Se copia la parte real del sistema de partida *) 


A:=Y; 
B:=G; (* Términos independientes *) 


(*-=Se incorporan a A los efectos de los diodos *) 
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For i:=1 to ND do 
with Di] do 
begin 
F := Id (Va); 
(* Corriente en esta iteración *) 
DF := Did (Va); 
(* Derivada de la corriente *) 


A[Na,Na] := A[Na,Na] + DE; 
(* Resistencia Equivalente *) 
A[Nc,Nc] := A[Nc,Nc] + DE; 
A[Na,Nc] A[Na,Nc] - DE; 
A[Nc,_Na] A[Nc,Nal] - DE; 


B[Nc] := B[Nc]+E - DE * Va; 

(* Generador de corriente equivalente *) 
BíNa] := Bí[Na]-E + DE * Va; 

end; 


(* Resolución del sistema en la presente iteración *) 
Resuelve _sistema(A,B,NN,Vi); 


(*- Se presentan los resultados *) 


Writeln; 
Writeln('Solución en la iteración ',NIt,':'); 
Presenta Solucion (Vi,NN); 


(* Se actualizan las tensiones en diodos y se 
comprueba si se cumple el márgen de error *) 


Final:=True; 
For i:= 1 to ND do 
with D[i] do 
begin 
F:=Vi[Na]-Vi[Nc]; 
(* Nuevo valor tensión en diodo i *) 
If Abs(Va-F) > E 
then Final:=False; 
(* Error diodo i mayor del permitido *) 


Va:= E; 
end; 
NIt:=NIt + 1; 
end; (* Del WHILE. Fin de las iteraciones *) 


(* Análisis de las causas de fin de las iteraciones *) 
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If Final 

then 
Writeln('El sistema converge”) 

else 
Nriteln('Tras las iteraciones fijadas el', 

"sistema no converge'); 
end; 
(* Del ELSE. Fin de tratamiento análisis no lineal *) 


end; 


(* Variables globales *) 
Var 
(* Matriz de admitancias del circuito *) 
Yn : Sistema; 
(* Matriz de generadores equivalentes *) 
leg : Vector; 
(* Número de nodos del circuito *) 
NNodos : Integer; 
(* Variables auxiliares *) 
lod : Integer; 
(* Matriz con la información de los diodos *) 
Matriz Diodos : Diodos; 
(* Número de diodos del circuito *) 
Numero Diodos : Integer; 
(* Error máximo permisible en el resultado final *) 


Error : Real; 
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(* Número máximo de iteraciones a realizar *) 


Max _Iteraciones : Integer; 


(* Procedimiento para visualizar la matriz de 
admitancias *) 


Procedure Ver_Matriz Admitancias(var Y: Sistema); 


var i,j: integer; 


begin 
Writeln; 
Writeln('Matriz de admitancias del circuito' == 
MA AA 5 
Writeln; 


for i:=1 to NNodos do 
for j:=1 to NNodos do 
AT A AA ASE 
end; 


(* Procedimiento para visualizar la matriz de 
generadores independientes *) 


Procedure Ver Generadores( var 1: Vector); 
var j:integer; 


begin 
Writeln; 
Writeln('Matriz de generadores independientes'); 
A == 
Writeln; 
for j:=1 to NNodos do 
AA E E 0 
end; 


(*= FIN DE LAS RUTINAS AUXILIARES DE VISUALIZACION = 
( A A 


(* Programa Principal *) 


BEGIN 


(* Borrado de la matriz de admitancias *) 
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for i:=0 to Max_Nodos do 
for j:=0 to Max_Nodos do 
Yn[i,j]:=0; 


(* Borrado matriz de generadores independientes *) 


for i:=0 to Max_Nodos do 
leq[i]:=0; 


(* Petición del número de nodos, comprobando si 
rebasa el máximo permitido *) 


repeat 
WriteLn('Introduzca el número de nodos del'); 
Write(' circuito (Sin contar el de masa)...>'); 


ReadLn(NNodos); 
If NNodos > Max_Nodos 
then writeln ('No es posible procesar tantos', 
' nodos con este programa! ); 
until NNodos <= Max_Nodos; 


Writeln; 


Writeln('ATENCION: Ponga su teclado en mayúsculas'); 
Writeln; 


(* Entrada de los componentes lineales del circuito *) 


Writeln('Introduzca los elementos lineales:'); 
Writeln; 


Entrada_del Circuito(Yn,leq); 
Ver_Matriz Admitancias (Yn); 


Ver_Generadores (leg); 


(* Entrada de los diodos *) 


WriteLn; 

repeat 
Write('Número de diodos en el circuito...... => 
ReadLn(Numero Diodos); 

until Numero _Diodos <= Max_Diodos; 


If Numero Diodos <> O 
then begin 
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A A AAA PAK A A A A PP —————— e 
A mero maximo dae 1TeraciOMes. ........“ )) 


A 
ReadiLniMax Iteraciones) 


A A A A 


Ejemplos 


Para ilustrar el uso de este programa se van a resolver dos circuitos. 
El primero es un simple circuito resistivo, conocido como atenuador 
en T puenteada. 


50 100 2 100 


15 50 


Fig. 35. Atenuador en T puenteada. 


A continuación se puede ver el proceso completo de ejecución del pro- 
grama con los resultados obtenidos: 


Introduzca el número de nodos del 
circuito (Sin contar el de masa)...>3 


ATENCION: Ponga su teclado en mayúsculas 
Introduzca los elementos lineales: 


Tipos posibles de componentes: 
- R: Resistencia. 
- V: Generador independiente de tensión. 
- Il: Generador independiente de corriente. 
-= Q: Fin del circuito 


Introduzca el tipo de componente..... >v 
Datos del generador: 

Terminal negativo conectado al nudo..>0 
Terminal positivo conectado al nudo..>1 


Valor - del Generado Es: saunas anar e > ES 
Valor del componente en serie........ >50 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>1 
Terminal negativo conectado al nodo..>2 


Valor del. COMpORENte. s.m. .nnaros e >100 
Controla algún generador?............ >N 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>2 
Terminal negativo conectado al nodo..>3 


Valor del ¡COMPORENTE. co. casa ciao >100 
Controla algún generador? ..o..o..oms.-. >N 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>2 
Terminal negativo conectado al nodo..>0 


Valor del Componenté. messi >25 
Controla algún generador? «ocioso... >N 
Introduzca el tipo de componente..... >R 


Datos del componente: 

Terminal positivo conectado al nodo..>1 
Terminal negativo conectado al nodo..>3 
valor del Componente. sisiseós eins me >10 


TZ 


Introduzca el tipo de componente 


Datos del componente: 


Terminal positivo conectado al nodo... 


Y(1,1)= 1 
Y(1,2) = -1 
MA) 2 
v(2,1)= -1 
Y(2,2)= 6 
Y(2,3) = -1 
Y(3,1) = -1 
3,3) = 1 
Y(3,3)= 1 


Matriz de generadores independientes 


e e e A 


Número de diodos en el circuito 


.3000000000E-01 
. 0000000000E-02 
. 0O000000000E-01 
. 0000000000E-02 
. 0000000000E-02 
. O0O00000000E-02 
. 0000000000E-01 
. 0000000000E-02 
.3000000000E-01 


1 = 3.0000000000E-01 
12) = 0.0000000000E+00 
= 0.0000000000E+00 


Solución del circuito: 

vl1 = 6.2771739130E+00 
vV2 = 1.8750000000E+00 
V3 = 4.9728260869E+00 


El siguiente ejemplo corresponde a un circuito no lineal con dos dio- 
dos como el de la figura siguiente: 


2 


Fig. 36. 


TS 


La ejecución del programa es como sigue: 


Introduzca el número de nodos del 
circuito (Sin contar el de masa)...>3 


ATENCION: Ponga su teclado en mayúsculas 
Introduzca los elementos lineales: 


Tipos posibles de componentes: 
- R: Resistencia. 
- V: Generador independiente de tensión. 
- I: Generador independiente de corriente. 
-= Q: Fin del circuito 


Introduzca el tipo de componente..... >V 
Datos del generador: 

Terminal negativo conectado al nudo.:>0 
Terminal positivo conectado al nudo..>1 


Valor del [JOneradobs 2. aiisja sión osas >3 
Valor del componente en serie........ >2 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>2 
Terminal negativo conectado al nodo..>0 


Valor del ¡COMPOnentesi veces ae 2) 
Controla algún generador? 0... a >N 
Introduzca el tipo de componente..... >R 


Datos del componente: 
Terminal positivo conectado al nodo..>3 
Terminal negativo conectado al nodo..>0 


Valor del COMPONEBNES-.; 5 oicaianao areas e >4 
Controla algún generador?............ >N 
Introduzca el tipo de componente..... >Q 


Y(L,1) = 5 .0000000000E-01 
MLSZ) = 0.0000000000E+00 
Y (1,3) = 0.0000000000E+00 
Y(2,1) = 0.0000000000E+00 
Y(2,2) = 3.3333333333E-01 
Y(2,3) = O.0000000000E+00 
Y(3,1) = 0.0000000000E+00 
Y(3,2) = 0.0000000000E+00 
Y(3,3) = 2.5000000000E-01 
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Matriz de generadores independientes 


1(1) = 1.5000000000E+00 
E(2):= 0.0000000000E+00 
1(3) = 0. 0000000000E+00 


Número de diodos en el circuito...... >2 
Diodo 1: 

Nodo: del ¡ANO O... cupido dd al ia > E 
Nodo: del ¡cátodo caras amas ea >2 
Tensión primera aproximación......... >0 
Diodo 2: 

Nodo! del ¡ANO 0. ae iia ala is >1 
Nodo del :cCAtOdO. ¿somera iaa a >3 
Tensión primera aproximación......... >0.5 
EXFEOFr MÁXIMO: PSrmisSliDle. estero mms >0.001 
Número máximo de iteraciones......... >20 


Solución en la iteración 1: 


v1 = 1.6499377684E+00 
v2 = 1.1550974902E+00 
V3 = 1.1599944758E+00 


Solución en la iteración 2: 


v1 = 1.6491434466E+00 
v2 = 1.1549186568E+00 
V3 = 1.1618215642E+00 
Solución en la iteración 3: 
v1 = 1.6491074077E+00 
v2 = 1.1548909201E+00 
V3 = 1.1619306242E+00 


El sistema converge 


En este caso particular la convergencia es muy rápida. Si se hubieran 
tomado otras condiciones iniciales, el sistema podía haber resultado diver- 


gente o bien oscilar en torno a unos valores arbitrarios. 


TO 


ANALISIS DE CIRCUITOS EN 
ALTERNA 


L análisis en alterna es el más sencillo de los tres que com- 
ponen el proceso completo. En él empleamos las impe- 
dancias de los componentes para obtener sus relaciones 
corriente-tensión, como se vio en el capítulo «Impedan- 
cias», en la página 20. 


El proceso que debemos seguir es el siguiente: 


1. Sustituir los componentes no lineales que existan 
por sus circuitos equivalentes para pequeña señal. Para 
===== ello deberemos conocer normalmente los resultados del 
análisis en continua, ya que algunos de los parámetros de este circuito de- 
penden de aquéllos. 

Con esto habremos llegado a un circuito con componentes exclusiva- 
mente lineales. 

2. Plantear las ecuaciones del circuito por el método del análisis no- 
dal. 

3. Resolver el sistema formado. 


El primer apartado se va a estudiar en detalle en este capítulo. El se- 
gundo ya se ha desarrollado con anterioridad, y para el tercero se puede 
consultar el apéndice «Resolución de sistemas de ecuaciones», en la pági- 
na 134, donde se trata uno de los posibles métodos de resolución junto 
con un par de ejemplos. Para un estudio más detallado sobre la resolución 
de sistemas de ecuaciones, el lector hará bien en consultar la bibliografía 
básica que encontrará al final del libro, ya que este tema se sale fuera de 
la extensión del presente. 

Del mismo modo que en el capítulo anterior, se incluye en éste un pro- 
grama para ilustrar los métodos matemáticos presentados. El programa es 
también completo, esto es, lleva incluidas todas las rutinas necesarias para 
su correcto funcionamiento. El lector que no esté interesado en la justifi- 
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cación matemática del procedimiento a seguir puede acudir directamente 
a dicho programa, ya que los comentarios que en él se incluyen pueden 
resultar casi suficientes para comprender su forma de trabajo. 


MODELO DE PEQUEÑA SEÑAL 
PARA EL TRANSISTOR BIPOLAR 


Llamamos a este modelo «de pequeña señal» porque, para su uso, es 
necesario que la señal alterna que se aplica al transistor no sea de gran am- 
plitud, para poder aproximar la respuesta del transistor en esa pequeña 
zona por una línea recta. 


Fig. 37. Modelo del transistor bipolar en pequeña señal. 


Normalmente se trabaja con transistores polarizados en zona activa, 
esto es, con la unión base-emisor en polarización directa y la unión base- 
colector en polarización inversa. El modelo de pequeña señal exige que la 
exploración de la señal alterna que se superpone a la continua sea tal que 
no lleve nunca el transistor al corte (ambas uniones polarizadas en inver- 
so) ni a la saturación (ambas uniones bajo polarización directa). 

El circuito al que nos referimos se halla en la figura 37. En él se han 
despreciado algunos efectos de poca importancia para resaltar los funda- 
mentales. Los parámetros que necesitaremos conocer son: 


e hfe: ganancia en corriente. Aunque no aparece en el circuito, la em- 
plearemos en el cálculo de otros componentes. Su valor viene dado por: 


A 


Ib 


LL 


Siendo Ic e Ib las del punto de polarización. 
e rbe: resistencia base-emisor. Viene dada por la expresión: 


rbe = =— > Vt = en = 0,025 a temperatura ambiente. 
q 


e Cbe, Cbc: Capacidades base-emisor y base-colector. Son de un valor 
muy pequeño (algunos picofaradios) y se pueden tomar de las caracterís- 
ticas del transistor dadas por el fabricante sin cometer un gran error. 

e rbc: resistencia base-colector. Se obtiene también de las caracterís- 
ticas dadas por el fabricante, considerando que 


rbe 


hre 


rbc = , hre: Ganancia de corriente inversa en emisor común. 


Su valor es muy alto (4 M2), por lo que normalmente se elimina. 
e gm: Transconductancia de salida. Su valor viene dado por: 
hfe 


e rce: resistencia de salida. Se obtiene a partir de las características del 
fabricante, ya que rce = 1/hoe. 


MODELO DE PEQUEÑA SEÑAL DE LOS 
TRANSISTORES DE EFECTO DE CAMPO 


Las ecuaciones que representan el comportamiento de los transistores 
FET de unión y de estructura MOS (MOSFET), en régimen de pequeña se- 
ñal, son del mismo tipo. Por tanto, los podremos sustituir por el mismo cir- 
cuito equivalente para esta fase de cálculo. 

La figura 38 presenta dicho circuito. 


Fig. 38. Modelo de los transistores de efecto de campo. 
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Sus parámetros se pueden derivar a partir de las hojas de característi- 
cas del fabricante. La variedad de notaciones con que se les denomina es 
muy grande. Aquí se presentan algunas de las más corrientes: 

gm, gfs, yfs: Transconductancia directa. 

gm0: Valor de gn cuando Vgs = 0. 

Vp, Vgs(ott): Tensión de estrangulamiento o de «Pinch-off». Es el valor 


de la tensión puerta-fuente por encima del cual la corriente drenador-fuen- 
te no aumenta. 


Idss: Corriente de drenador para Vgs = 0. 


Ciss, Cgss: Capacidad de entrada en configuración de fuente común, 
con la salida en cortocircuito. 


Crss: Capacidad inversa de transferencia, con la entrada cortocircuitada. 


Coss: Capacidad de salida en fuente común, con la entrada cortocircui- 
tada. 


gd, yos, gos: Conductancia de salida en fuente común. 


Las ecuaciones que podemos aplicar son: 
gm = gm0 (1-85, 
Vp 
2 Idss 


Vp = 
E gm0 


Cds = Coss — Crss Normalmente es muy pequeña (<1pf) y se desprecia. 
1 
-— V-Vgd 
ue 
V-Vgs 


Cgd = Crss o bien, aproximadamente, Cgd 
Cgs = Ciss — Crss o bien, aproximadamente, Cgs = 


Es preciso señalar, sin embargo, que los datos del fabricante suelen te- 
ner una dispersión tal que hacen que el modelo sea solamente aproximado. 


MODELO DEL AMPLIFICADOR OPERACIONAL 


El amplificador operacional es un circuito integrado lineal, común- 
mente empleado hoy en día como un elemento más de circuito, sin ocu- 
parse de su estructura interna. Es por esto por lo que resulta interesante 
disponer de un modelo de él, que nos permita realizar el análisis de los cir- 
cuitos en los que intervenga de una forma sencilla. 
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La representación simbólica del amplificador operacional se ofrece en 
la figura 39. 


v+ 


vV- 
Fig. 39. Diagrama del amplificador operacional. 


Cuando se analizan circuitos con amplificadores operacionales por los 
métodos tradicionales (es decir, con lápiz y papel), se realizan una serie 
de simplificaciones para facilitar el cálculo. En el análisis por ordenador 
no es necesario realizarlas, con lo que conseguimos un estudio más próxi- 
mo a la realidad. En el siguiente cuadro se van a estudiar las principales 
diferencias entre el análisis simplificado y un modelo más real: 


_ Parámetro | Simb. | Mod. simplif. | Mod. real 


Con estas consideraciones podemos establecer el siguiente modelo del 
amplificador operacional para el análisis en el dominio de la frecuencia: 


Fig. 40. Modelo de un amplificador operacional real en pequeña señal. 


Los parámetros Av, Rin y Ro se obtienen directamente de las caracte- 
rísticas dadas por el fabricante. En cuanto a R y C se fijan a partir de la 
frecuencia dada por el fabricante para el polo dominante en circuito abier- 
to, que denominaremos f0, según la relación: 


1 


R.C = ———— 
DATO 
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Valores aproximados de estos parámetros son: 


Ri =2MQ Av = 200.000 
Ro =75 Q fO = 10 Hz. 


Con esto hemos logrado simplificar un circuito complejo, que ocupa- 
ría un volumen de memoria en nuestro ordenador enorme y un tiempo de 
cálculo bastante extenso, por otro más simple y que, con pequeñas restric- 
ciones, se comporta prácticamente igual para señales alternas de baja am- 
plitud. Esta técnica recibe el nombre de macromodelado, y es extensible 
a Cualquier otro tipo de circuito de uso general. 

El macromodelo que realicemos puede ser tanto más general como 
nuestro análisis lo requiera. En cualquier caso, deberá reflejar los princi- 
pales fenómenos que ocurran en la realidad bajo las condiciones de tra- 
bajo. En el caso del amplificador operacional, por ejemplo: 


e Para el análisis en continua, se incluirían: 


Corrientes de polarización. 

Tensiones de offset. 

Característica de transferencia estática. 
Derivas de los parámetros. 


e Para el análisis de pequeña señal más detallado: 


e Capacidad de entrada. 
e Factor de rechazo del modo común (CMMR). 
e Figura de ruido. 


e Para el análisis en gran señal: 
e Slew-rate (velocidad de subida). 


e Para el análisis transitorio, en función de la respuesta a un escalón 
de tensión en su entrada: 


e Tiempo de subida. 
e Tiempo de asentamiento. 
e Tiempo de adquisición. 


Que son de interés en muchas aplicaciones (comparadores, muestreo-re- 
tención). 


EJEMPLO DE PROGRAMACION 


Continuando en la línea de los ejemplos hasta ahora introducidos, se 
ofrece este nuevo programa. En él se emplean las rutinas de entrada del 
capítulo «Fundamentos del diseño de circuitos asistido por ordenador», en 
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la página 13, del mismo modo que las rutinas de resolución de un sistema 
de ecuaciones complejas presentadas en el «Resolución de sistemas de 
ecuaciones», en la página 134. 

Como se puede apreciar, el programa no tiene la complejidad del caso 
de resolución de circuitos no lineales en continua. Por el contrario, los re- 
sultados que de él se obtienen no dependen de parámetros a fijar de forma 
aproximada «a priori», con lo que su fiabilidad es mucho mayor. 

De igual forma que en anteriores programas, los frecuentes comenta- 
rios deben ser suficientes para comprender el funcionamiento de las ruti- 
nas que aquí nos ocupan. En este caso, su brevedad hace innecesaria cual- 
quier otra posible explicación. Pasemos, pues, a verlo: 


A A A A 
¡AX NOdCdOSj Or reaj 
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begin 
writeln('Datos del generador:'); 
write('Terminal negativo conectado al nudo..>'); 


readln(N1); 

write('Terminal positivo conectado al nudo..>'); 
readln(N2); 

E AAA => 
readln(valor); 

write('Valor del componente en serie........ >; 


readlin(G serie); 
G_Serie := 1/G_Serie; 
Valor:=Valor * G_Serie; 
end; 
procedure Gen Corriente(var N1,N2: integer; 
var valor : real); 


begin 
writeln('Datos del generador:'); 
write('Terminal origen conectado al nudo....>'); 
readln(N1); 


write('Terminal destino conectado al nudo...>'); 
readln(N2); 
write('ValorProgram Analisis_en_alterna; 


(* Este programa realiza el análisis en alterna de 
un circuito para un margen de frecuencias 
especificado *) 


(*=====================o=o==========================* 
(*= PROCEDIMIENTOS PARA LA ENTRADA DEL CIRCUITO =*) 
=== _ _z-- o e e _———————— * ) 


(* Estos procedimientos son los mismos que los 
correspondientes del programa Entrada_de Circuito, 
ya estudiado anteriormente, y por eso se presentan 
aquí sin comentarios. *) 

Const 

Max_Nodos = 20; 
Type 

Circuito  —= Array [0..Max_Nodos,0..Max_Nodos, 

: =11)]-of real; 

Generadores = Array[0..Max_Nodos] of real; 

Procedure Entrada del Circuito ( var Y: Circuito; 
var Il: Generadores); 

var 

Componente : Char; 

N1, 
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TipoG : Integer; 


ValorComponente, 
ValorGenerador, 

GGenerador : Real; 
Error : Boolean; 


procedure Gen Tension(var N1,N2: integer; 
var valor,G_ serie : real); 


begin 


writeln('Datos del generador:'); 
write('Terminal negativo conectado al nudo... 
readln(N1); 

write('Terminal positivo conectado al nudo.. 
readln(N2); 

A A 
readin(valor); 

write('Valor del componente en serie........ 
readln(G_serie); 

G-Serie :=-1/6_Serie; 

Valor:=Valor * G Serie; 


end; 
procedure Gen Corriente(var Nl1,N2: integer; 


var valor : real); 


begin 


writeln('Datos del generador:'); 
write('Terminal origen conectado al nudo... 
readln(N1); 

write('Terminal destino conectado al nudo... 
readln(N2); 

A 
readin(valor); 


end; 
procedure Entrada Componente (C: integer); 
var 


Control : Char; 


begin 
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TipoG:=0; 

Error:=False; 

writeln('Datos del componente:'); 
write('Terminal positivo conectado al nodo.. 
readin(N1); 

write('Terminal negativo conectado al nodo. 
readln(N2); 

write('Valor del componente................. 
readin(ValorComponente); 


if ((c= 0) or (c = -1)) then ValorComponente:= 
1/ValorComponente; 


Control:='" * 


while not === in-EN5,%5851)Y 40 
begin 


readlin(Control); 


end; 
if Control ='S' then 
begin 
write('Variable que lo controla (1/V)....... == 


readln(Control); 

case control of 
'v': TipoG:=1; 
“: TipoG:=2; 
else Error:=True; 


end; 
write('Tipo de generador (1I/V).............. su; 
readin(Control); 
case Control of 
"v!': begin 
TipoG:=TipoG + 10; 
Gen Tension(NC1, NC2, ValorGenerador, 
GGenerador); 
end; 
'I1': begin 
TipoG:=TipoG + 20; 
Gen Corriente(NC1, NC2,ValorGenerador); 
end; 
else 
Error:=True; 
end; 
end; 


If not Error then 
begin 
Y[N1,N1,C]:=Y[N1,N1,C]+ValorComponente; 
Y[N2,N2,C]:=Y[N2,N2,C]+ValorComponente; 
Y[N1,N2,C]:=Y[N1,N2,C]-ValorComponente; 
Y[N2,N1,C]:=Y[N2,N1,C]-ValorComponente; 
if TipoG div 10 = 1 then 
begin 
Y[NC1,NC1,0]:=Y[NC1,NC1,0]+GGenerador; 
Y[NC2,NC2,0]:=Y[NC2,NC2,0]+GGenerador; 
Y[NC1,NC2,0]:=Y[NC1,NC2,0]=-=GGenerador; 
YÍ[NC2,NC1,0]:=Y[NC2,NC1,0]-GGenerador; 
end; 
if TipoG mod 10 = 1 then 
begin 
Y[NC1,N1,0]:=Y[NC1,N1,0]+ValorGenerador; 
Y[NC2,N2,0]:=Y[NC2,N2,0]+ValorGenerador; 
Y[NC1,N2,0]:=Y[NC1,N2,0]-ValorGenerador; 
Y[NC2,N1,0]:=Y[NC2,N1,0]-ValorGenerador; 
end; 
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if TipoG mod 10 = 2 then 
begin 
Y[NC1,N1,0]:=Y[NC1,N1,0]+ValorGenerador* 


ValorComponente; 


Y[NC2,N2,0]:=Y[NC2,N2,0]+ValorGenerador* 


ValorComponente; 


Y[NC1,N2,0]:=Y[NC1,N2,0]-ValorGenerador* 


ValorComponente; 


Y[NC2,N1,0]:=Y[NC2,N1,0]-ValorGenerador* 


end; 
end 


ValorComponente; 


else writeln('Entrada equivocada. 


end; 
begin 
writeln('Tipos posibles de componentes:'); 
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! 
'Repítala, por favor.'); 


writeln(' - R: Resistencia.'); 
writeln(' - L: Bobina.'); 
writeln(' - C: Condensador.'); 
writeln(' - V: Generador indep. de tensión. == 
writeln(' - 1: Generador indep. de corriente.'); 
writeln(' - Q: Fin del circuito' ); 
repeat 
WriteLn; 
Write('Introduzca el tipo de componente..... E 


ReadEn(Componente); 
Case Componente of 


a 
A 
ra 
q! 


end; 


Entrada_Componente(0); 


: Entrada Componente(-1); 
: Entrada Componente(1); 
: begin 


Gen Tension(NC1, NC2, ValorGenerador, 
GGenerador); 
Y[NC1,NC1,0]:=Y[NC1,NC1,0]+GGenerador; 
Y[NC2,NC2,0]:=Y[NC2,NC2,0]+GGenerador; 
Y[NC1,NC2,0]:=Y[NC1,NC2,0]=-GGenerador; 
Y[NC2,NC1,0]:=Y[NC2,NC1,0]-GGenerador; 
I[NC1]:=I[NC1]-ValorGenerador; 
I[NC2]:=I[NC2]+ValorGenerador; 
end; 


: begin 


Gen Corriente(NC1, NC2, ValorGenerador); 
I[NC1]:=1[NC1]- -ValorGenerador; 
I[NC2]:=I[NC2]+ValorGenerador; 

end; 


until Componente='0Q' 


end; 


(* Estos procedimientos son los mismos que los 
correspondientes del programa 
Sistema de Ecuaciones Complejas, presentado en los 
apéndices, y por eso se presentan aquí sin 
comentarios. *) 


Const 


Orden Maximo = 20; 


Type 
Complejo= 
Sistema = 


Vector = 
Procedure R 
var 

R: real; 
begin 


array[1..2] of real; 
array[1..Orden_Maximo,1..Orden- Maximo] 

of Complejo; 
array[1..Orden _ Maximo] of complejo; 
P(Var C:Complejo):; 


R:=Sqgrt(Sgr(C[1])+Sgr(C12])); 
if Abs(C[1]) < 1E-10 


then 
begin 
if Ccr2]>0 
then C[2]:=Pi/2 
else C[2]:=-Pi/2; 
end 
else 
begin 
CERO 
then C[2]:=ArcTan(C[2]/C11)) 
else C[2]:=ArcTan(Cc[2]/C[1]) + Pi; 
end; 
CEHE=R> 
end; 


Procedure Iguc(C1: Complejo;var Cr:Complejo); 


var 


i:integer; 


Begin 


For -i:=1-to-2 do Crtik:=C1iE]3; 


end; 
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Procedure SumC(C1,C2: Complejo;var Cr:Complejo); 
var 
i:integer; 
Begin 
For i:=1 to 2 do Cr[i]:=C1[i]+C2[i]; 
end; 
Procedure ResC(C1,€2: Complejo;var €r:Complejo); 
var 
i:integer; 
Begin 
For i:=1 to 2 do CrliJ:=C1ftiJ-C2ti]H; 
end; 
Procedure MuiC(C1,C€2: Complejo;var-C€r:Complejo); 
var 
ivinteger; 
Begin 
AEREA EEE E ao 
CcrR2j=CH3*C212]4+C112]3*C0211]; 
end; 
Procedure DivC(C1,C2: Complejo;var Cr:Complejo); 
var 
izinteger; 
D: real; 
Begin 
D:=sqr(C2[1])+Sqr(Cc2[2]); 
Crf1]:=(C1[11]1*C2[1]+C112]*C212])/D; 
ert2]:=(-C111]*C0212]+0112]*0211)7D; 
end; 
Procedure Resuelve Sistema(var A: Sistema; 
var B: vector; Orden: Integer;var X:Vector); 


var 
1,3,k-— integer; 
>= Vector; 
[el : Complejo; 
Begin 
For-—j:=1-to-Orden- do 
Begin 
for i:=j to Orden do 
begin 
for k:=1-to-j-1-do 
begin 
MulC(A[i,k],A[k,3],C); 
ResC(A[i,3j)],C,A[i,31); 
end; 
end; 
fori:=j*1- to Orden do 
begin 


for k:=1 to j-1 do 
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begin 
MICA ALA iO; 
ResC(A(5,1],C,46,1)); 
end; 
DIVA EALABGAB iD 
end; 
end; 


For-i:=1 to Orden do 
begin 
IgucBhYA LED; 
for—k:=1-to-i-3-do 
begin 
MulC(Afixk],YEx3,C); 
ResC(Y[i],C,Y[i]); 
end; 
DICOAITAEAEAAADO 
end; 
For—i:=0rden downto-1-do 
begin 
ge 
for—k:=0rden-—downto-3+1-do 
begin 
MulC(A[i,k],X[k],C); 
ResC(X[i],C,X[fi1); 
end; 


(*====F TN RESOLUCION-—DE-—UN—SISTEMA—DE—ECUACIONES====X)> 


(* Procedimiento para resolver el circuito en alterna 
a una frecuencia determinada. Parámetros de 


entrada: 
=- Y ¡Matriz de admitancias del-circuito: 
-1 : Matriz de generadores de corriente. 
=N : Numero de nodos, 
=P Frecuencia a la que se efectúa el análisis, 
=Y Solución del- circuito a la frecuencia dada 


(vector complejo).*) 


Procedure Solucion-en F(VarY:—circuito; 
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Var Ig:Generadores; N: integer; 
F: real; Var V:Vector); 


Var 


Ya : Sistema; 

(* Matriz del sistema de ecuaciones a una frecuencia*) 
B-—Vector; (* Términos independientes *) 
i,j: integer; 


Begin 
E:=2*Pi*E; 
(* Cálculo de la pulsación w a partir de la frec. *) 


(*A-continuación se calcula la matriz de admitancias 
del sistema para la frecuencia de interés *) 


For i:=1 to N do 


begin 
For-3:=1-to N do 
begin 
Yat1,3H1]:=Y141, 3,01]; 
AAA AAA 
end; 
Bfi111]:=Ig(fi); 
AROS 
end; 


(* Se resuelve el sistema obtenido *) 


Resuelve Sistema(Ya,B,N,V); 
end; 


(* Variables globales *) 


Var 
Yn ECTS 
(* Matriz de admitancias del circuito *) 
leg : Generadores; 


(* Matriz de generadores equivalentes *) 

NNodos “Integer; 
(* Número de nodos del circuito *) 
E - (*_ Frecuencia a analizar *) 
Fmax, (* Frecuencia máxima a analizar *) 
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dE : Real; (* Incremento de frecuencia *) 
Vk : Vector; (*=Solución del circuito *) 
LE : Integer; (* Variables auxiliares *) 


Procedure Ver Matriz Admitancias(var Y: Circuito); 


var i1,j: integer; 


begin 
Writeln; 
Nriteln('*Matriz de admitancias del circuito” == 
MA AAA = Y 
Writeln; 


for i:=1 to NNodos do 
for j:=1- to NNodos do begin 


A A A A A A A A 
Witein?-—=== 1H 4955535835 
end; 


end; 

Procedure Ver Generadores( var I: Generadores); 
var j:integer; 

begin 


Writeln; 
Writeln('Matriz de generadores independientes'); 
Mide =7 
Writeln; 
for 3j:=1- to NNodos- do 
WriteLn('1(*,3,') = *,113)); 
end; 


AE O AE RR O 


(* Programa Principal *) 
BEGIN 


(* Borrado de la matriz de admitancias *) 
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for-1:=0-to Max Nodos-do 
for—3:=0-to Max _Nodos-do 
for-k:==1-to-1-do 
Yn[i,j,kx]:=0; 


(*—Borrado matriz de generadores independientes *) 


for i:=0 to Max_Nodos do 
legq[i]:=0; 


(Petición del número de nodos, comprobando si 


rebasa el máximo permitido *) 


repeat 
Writebn(' Introduzca el número de nodos del'); 
Nritef'—circuito (Sin—contar-el de masa)...>'); 


ReadLn(NNodos); 
If NNodos > Max_Nodos 
then writeln ('No es posible procesar tantos', 
'—nodos con este programa”); 

until-NNodos-<=-Max-Nodos; 
Writeln; 
Writeln('ATENCION: Ponga su teclado en mayúsculas”); 
WNriteln; 


Entrada del Circuito(Yn,leg); 
Ver_Matriz_Admitancias(Yn); 


Ver _Ceneradores(leg); 


WriteLn; 
WriteLn(' ANALISIS EN ALTERNA DEL CIRCUITO" == 
A ei == 
WriteLn; 


(* Se piden los límites de frecuencia para los que se 
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debe realizar el análisis *) 


Write("Erecuencia mínima para analizar...... >); 
ReadEntE)3; - 

Write('Frecuencia máxima para analizar...... E 
ReadLn(Emax); 

Write(*Factorincremento de frecuencias... 3 
ReadEn(dF); 


While F-<=-Emax-do 
begin 


Para ilustrar su funcionamiento se va a analizar una etapa amplificado- 
ra en emisor común con acoplo en alterna, cuyo esquema eléctrico es 
como sigue: 


Fig. 41. Etapa amplificadora en emisor común. 
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El circuito equivalente, empleando el modelo del transistor para pe- 
queña señal, tiene la forma: 


nO Cc, 


g 
Fig. 42. Circuito equivalente en pequeña señal. 


En el que se ha puesto ya la asignación de nodos elegida para realizar el 
análisis. 
Los valores de los componentes son: 


Va =1 RE =1002 Cl =1pF 
Rg =50Q rbe = 1KQ  C2=1puF 
RB =RBl //RB2=25KQ rce = 35KQ gm =0,36 
Re =1 KQ Cbc = 10 pF 
R1 =1KQ Cbe = 25 pF 


Los valores de los parámetros del transistor se obtienen de las condi- 
ciones fijadas por su punto de polarización, como se ha explicado en este 
capítulo. 

Proceso y resultado del análisis vienen resumidos en el siguiente lis- 
tado: 


Introduzca el número de nodos del 
circuito (Sin contar el de masa)...>5 
ATENCION: Ponga su teclado en mayúsculas 


Tipos posibles de componentes: 


- R: Resistencia. 
- L: Bobina. 
- C: Condensador. 
- V: Generador indep. de tensión. 
- Il: Generador indep. de corriente. 
-= Q: Fin del circuito 
Introduzca el tipo de componente..... >V 


Datos del generador: 

Terminal negativo conectado al nudo..>0 
Terminal positivo conectado al nudo..>1 
Valor del gene bado. 0. acta e >1 
Valor del componente en serie........ >50 
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Introduzca el tipo de componente..... 
Datos del componente: 

Terminal positivo conectado al nodo.. 
Terminal negativo conectado al nodo... 
Valor del. COMPOnentSe.. mimi easier 


Controla algún generador?............ 


Introduzca el tipo de componente..... 
Datos del componente: 

Terminal positivo conectado al nodo.. 
Terminal negativo conectado al nodo... 
valor del componente......c0.0... <<... 
Controla algún generador?...........-. 


Introduzca el tipo de componente..... 
Datos del componente: 

Terminal positivo conectado al nodo.. 
Terminal negativo conectado al nodo.. 


Variable que lo controla (I/V)....... 
Tipo de generador (I/V)............-.- 
Datos del generador: 

Terminal origen conectado al nudo.... 
Terminal destino conectado al nudo... 
valor dell GENESFAÍOE.. -. e sioicio sea a es 


Introduzca el tipo de componente..... 
Datos del componente: 

Terminal positivo conectado al nodo.. 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente..... 


Datos del componente: 
Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo... 


Introduzca el tipo de componente..... 
Datos del componente: 

Terminal positivo conectado al nodo.. 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente..... 
Datos del componente: 
Terminal positivo conectado al nodo.. 


>N 


>0.36 


>25E-12 
>N 
>C 


>2 
>4 
>10E-12 
>N 
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Terminal negativo conectado al nodo... 


Introduzca el tipo de componente.... 


Datos del componente: 


Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo... 


Introduzca el tipo de componente.... 


Datos del componente: 


Terminal positivo conectado al nodo. . 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente.... 


Datos del componente: 


Terminal positivo conectado al nodo... 
Terminal negativo conectado al nodo.. 


Introduzca el tipo de componente.... 


Matriz de admitancias del circuito 


Y(1,1) =  0.0000000000E+00 1/jw 


+ 2 .0000000000E-02 + 1.0000000000E-06 


Y(1,2) = 0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + -1.0000000000E-06 


Y(1,3) =  0.0000000000E+00 1/jw 


sd 0.0000000000E+00 + 0.0000000000E+00 


Y(1,4) =  0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + 0.0000000000E+00 


Y(1,5) = 0.0000000000E+00 1/jw 
+ 


0.0000000000E+00 + O.0000000000E+00 


Y(2,1) =  0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + -1.0000000000E-06 


LUZ 2) == 0.0000000000E+00 1/jw 


+ 1.0400000000E-03 + 1.0000350000E-06 


Y(2,3) =  0.0000000000E+00 1/jw 


+  -1.0000000000E-03 + -2.5000000000E-11 


Y(2,4) =  0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + -1.0000000000E-11 


Y(2,5) =  0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + 0. 0000000000E+00 


v(3,1) 


0.0000000000E+00 1/jw 


+ 0.0000000000E+00 + 0. 0000000000E+00 


jw 
jw 
jw 
jw 
jw 
jw 
jw 
jw 
jw 
jw 


jw 


Y(3,2) = 0.0000000000E+00 1/jw 

+  -3.6100000000E-01 + -2.5000000000E-11 jw 
Y(3,3) =  0.0000000000E+00 1/jw 

+  3.7102857143E-01 +  2.5000000000E-11 jw 
Y(3,4) =  0.0000000000E+00 1/jw 

+ -2.8571428571E-05 +  0.0000000000E+00 jw 
Y(3,5) =  0.0000000000E+00 1/jw 

+  0.0000000000E+00 +  0.0000000000E+00 jw 
Y(4,1) = 0.0000000000E+00 1/jw 

+  0.0000000000E+00 +  0.0000000000E+00 jw 
Y(4,2) = 0.0000000000E+00 1/jw 

+  3.6000000000E-01 + -1.0000000000E-11 jw 
Y(4,3) = 0.0000000000E+00 1/jw 

+ -3.6002857143E-01 +  0.0000000000E+00 jw 
Y(4,4) =  0.0000000000E+00 1/jw 

+  1.0285714286E-03 +  1.0000100000E-06 jw 
Y(4,5) = 0.0000000000E+00 1/jw 

+  0.0000000000E+00 + -1.0000000000E-06 jw 
Y(5,1) = 0.0000000000E+00 1/jw 

+  0.0000000000E+00 +  0.0000000000E+00 jw 
Y(5,2) = 0.0000000000E+00 1/jw 

+  0.0000000000E+00 +  0.0000000000E+00 jw 
Y(5,3) = 0.0000000000E+00 1/jw 

+  0.0000000000E+00 +  0.0000000000E+00 jw 
Y(5,4) =  0.0000000000E+00 1/jw 

+  0.0000000000E+00 + -1.0000000000E-06 jw 
Y(5,5) = 0.0000000000E+00 1/jw 

+  1.0000000000E-03 +  1.0000000000E-06 jw 


Matriz de generadores independientes 


1(1) =  2.0000000000E-02 
1(2) =  0.0000000000E+00 
1(3) =  0.0000000000E+00 
1(4) =  0.0000000000E+00 
1(5) =  0.0000000000E+00 


Frecuencia mínima para analizar...... >1 
Frecuencia máxima para analizar...... >10E6 
Factor incremento de frecuencias..... >10 


Solución en F=—  1.0000000000E+00 Hz. 
V1: 9.9997E-01 +j( -3.115E-04) 

==> 9.9997E-01 * explj(-1.785E-02)] 
V2: 8.5442E-03 +j( 9.1911E-02) 

==> 9.2308E-02 * explj(8.4689E+01)] 
V3: 8.3064E-03 +j( 8.9359E-02) 

==> 8.9744E-02 * exp[j(8.4689E+01)] 
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VA: -8.842E-02 +j( 
==> 8.9482E-01 * 
V5: 5.5911E-03 +j( 
==> 5.6222E-03 * 


Solución en F= a 
V1: 9.9843E-01 +j( 
==> 9.9843E-01 * 
V2: 4.6234E-01 +j( 
==> 6.7902E-01 * 
V3: 4.4947E-01 +j( 
==> 6.6016E-01 * 
V4: -4.745E+00 +j( 
==> 6.5443E+00 * 
VS: 2.6339E-01 +j( 
==> 4,1038E-01 * 


Solución en FE= A. 


Vl: 9.9667E-01 +j( 
==> 9.9667E-01 * 
V2: 9.8552E-01 +j( 
==> 9.9116E-01 * 
V3: 9.5835E-01 +j( 
==> 9.6386E-01 * 


V4: -6.881E+00 +3( 
==> 7.0680E+00 * 
V5: -2.676E+00 +j( 


==> 3.7603E+00 * 


Solución en F= al 
V1: 9.9664E-01 +j3( 
==> 9.9664E-01 * 
V2: 9.9647E-01 +j( 
==> 9.9653E-01 * 
V3: 9.6916E-01 +j( 
==> 9.6922E-01 * 


V4: -4.867E+00 +j( 
==> 4.8778E+00 * 
V5: -4.798E+00 +j( 
==> 4.8172E+00 * 
Solución en F= Le 
V1: 9.9664E-01 +j( 


==> 9.9664E-01 * 
V2: 9.9658E-01 +3( 
==> 9.9658E-01 * 


V3: 9.6927E-01 +j( 
==> 9.6927E-01 * 
VA: -4.833E+00 +j( 


==> 4.8332E+00 * 
VS: -4.832E+00 +3( 
==> 4.8325E+00 * 


-8.904E-01) 
exp[j(2.6433E+02)] 
-5.907E-04) 

exp[j(-6.031E+00)] 


. 0000000000E+01 Hz. 


-1.684E-03) 
exp[j(-9.665E-02)] 
4.9730E-01) 
exp[j(4.7087E+01)] 
4.8351E-01) 
exp[j(4.7089E+01)] 
-4.507E+00) 
exp[j(2.2352E+02)] 
-3.147E-01) 
exp[j(-5.007E+01)] 


0000000000E+02 Hz. 
-3.503E-04) 
exp[j(-2.014E-02)] 
1.0564E-01) 
exp[j(6.1185E+00)] 
1.0291E-01) 
exp[j(6.1292E+00)] 
1.6159E+00) 
exp[j(1.6678E+02)] 
-2.642E+00) 
exp[j(2.2464E+02)] 


. 0000000000E+03 Hz. 


-5.302E-05) 
exp[j(-3.048E-03)] 
1.0640E-02) 
exp[j(6.1177E-01)] 
1.0378E-02) 
exp[j(6.1351E-01)] 
3.3021E-01) 
exp[j(1.7612E+02)] 
-4.334E-01) 
exp[j(1.8516E+02)] 


0000000000E+04 Hz. 
-1.881E-04) 
explj(-1.081E-02)] 
8.8133E-04) 
exp[j(5.0670E-02)] 
8.6040E-04) 
exp[j(5.0860E-02)] 
3.6016E-02) 
exp[j(1.7957E+02)] 
-4.089E-02) 
exp[j(1.8048E+02)] 


Solución en F= 1.0000000000E+05 Hz. 
V1: 9.9663E-01 +j( -1.847E-03) 

==> 9.9663E-01 * explj(-1.062E-01)] 
V2: 9.9657E-01 +j( -1.739E-03) 

==> 9.9657E-01 * exp[lj(-1.000E-01)] 
V3: 9.6926E-01 +j( -1.689E-03) 

==> 9.6927E-01 * explj(-9.984E-02)] 
V4: -4.833E+00 +j( 3.0792E-02) 

==> 4.8326E+00 * exp[j(1.7963E+02)] 
VS: -4.833E+00 +j( 2.3101E-02) 

==> 4.8326E+00 * exp[j(1.7973E+02)] 


Solución en F= 1.0000000000E+06 Hz. 
Vl1: 9.9572E-01 +j( -1.842E-02) 

==> 9.9589E-01 * exp[j(-1.060E+00)] 
V2: 9.9566E-01 +j( -1.840E-02) 

==> 9.9583E-01 * exp[j(-1.059E+00)] 
V3: 9.6838E-01 +j( -1.787E-02) 

==> 9.6854E-01 * exp[j(-1.057E+00)] 
V4: -4.819E+00 +j( 2.7431E-01) 

==> 4.8268E+00 * explj(1.7674E+02)] 
V5: -4.819E+00 +j( 2.7354E-01) 

==> 4.8268E+00 * exp[j(1.7675E+02)] 


Solución en F= 1.0000000000E+07 Hz. 
V1: 9.2285E-01 +j( -1.478E-01) 

==> 9.3461E-01 * explj(-9.099E+00)|] 
V2: 9.2280E-01 +j( -1.478E-01) 

==> 9.3456E-01 * exp[j(-9.098E+00)] 
V3: 8.9759E-01 +j( -1.435E-01) 

==> 9.0899E-01 * exp[j(-9.083E+00)] 
V4: -3.735E+00 +j( 2.1985E+00) 

==> 4.3339E+00 * exp[j(1.4952E+02)] 
V5: -3.735E+00 +j( 2.1984E+00) 

==> 4.3339E+00 * exp[j(1.4952E+02)] 


La tensión de mayor interés es la del nodo 5, que es la salida del cir- 
cuito. Si se estudia su valor modular (tercera columna de los resultados 
de salida), se puede ver como hasta pasados los 100 Hz no alcanza su valor 
estable (del orden de 4.8 V.). Esto, como es natural, es debido al acoplo 
en alterna, realizado a través de los condensadores C1 y C2. 


La ganancia se puede calcular para cada frecuencia, realizando la divi- 
sión V5/V1. Vemos cómo es aproximadamente del orden de 5 en las fre- 
cuencias centrales, con una fase de unos 180 grados (es un amplificador 
inversor), dato obtenido de la cuarta columna de los resultados. 


Hasta muy alta frecuencia no se hacen notar los efectos de las capaci- 
dades internas del transistor Cbe y Cbc. Incluso a 10 MHz, máxima fre- 
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cuencia de nuestro análisis, la disminución de la salida es sólo de 0.5 V., 
mientras que en la fase ya es de 30 grados. 

Como puede verse, el análisis en frecuencia es quizá una de las partes 
más sencillas del proceso de análisis de un circuito. No por ello es la me- 
nos importante. Normalmente su importancia es decisiva, sobre todo en 
caso de amplificadores y otros circuitos en los que se busca una respuesta 
de tipo lineal. El conocimiento a fondo de sus métodos de análisis es, por 
tanto, fundamental. 
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ANALISIS DE CIRCUITOS EN 
REGIMEN TRANSITORIO 


N el análisis transitorio se busca calcular la evolución de 
la respuesta de un circuito, desde un instante de tiempo 
al siguiente, cuando a éste se le aplica una excitación re- 
pentinamente. Ñ 

En caso de que el circuito sea lineal, el funcionamien- 
to puede describirse por un conjunto de ecuaciones dife- 
renciales, que, si los parámetros no varían con el tiempo, 
serán lineales y con coeficientes constantes. 

Si los circuitos son no lineales, su descripción vendrá 
dada, en general, por un conjunto de ecuaciones diferenciales no lineales. 
En este caso los métodos de resolución analíticos son inútiles y se preci- 
san técnicas numéricas. 

Los métodos comúnmente empleados para atacar este problema se ba- 
san en el planteamiento de las ecuaciones de estado del circuito. Estas 
ecuaciones son las que relacionan entre sí las variables de estado de la red, 
que no son más que las corrientes en las bobinas y tensiones en los con- 
densadores. Se denominan así porque sus derivadas con respecto del tiem- 
po especifican cómo el «estado» en que se encuentra el circuito puede pre- 
decirse desde un instante de tiempo al siguiente. 

Una vez que se tienen estas ecuaciones de estado se resuelve el sistema 
de ecuaciones diferenciales por un método numérico, tanto más complejo 
cuanto más alta sea la exactitud que necesitamos en las soluciones. 

Otro procedimiento se basa en aplicar las fórmulas de integración nu- 
mérica a las ecuaciones diferenciales que definen condensadores y resis- 
tencias, esto es, a sus relaciones corriente-tensión. La solución de esta in- 
tegración numérica dará lugar a un circuito equivalente que podremos sus- 
tituir en la red, y proceder así al análisis en continua de un circuito muy 
simplificado. . 

Este es el método denominado del circuito incremental: Es el que se va 
a desarrollar en este capítulo, dejando los otros procedimientos para los 
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lectores de profundos conocimientos matemáticos que se sientan atraídos 
por la aventura del cálculo numérico. En la bibliografía citada al final de 
la obra encontrarán suficiente material para explorar en sus ratos de ocio. 


ANALISIS TRANSITORIO DE CIRCUITOS LINEALES 


A continuación vamos a realizar la integración de las ecuaciones 
corriente-tensión de los componentes básicos: condensadores, bobinas y 
resistencias. 

Podemos emplear cualquier método de integración numérica. Cuanto 
mayor sea la precisión requerida, más complejo será el método a utilizar. 
En todo el desarrollo que sigue se emplea el método del trapecio, que el 
lector profano puede encontrar explicado con detalle en el apéndice «Re- 
solución de una ecuación diferencial por el método del trapecio», en la pá- 
gina 131. 


1. Condensador 
Despejando la derivada en la ecuación del condensador se obtiene: 


A 
dt C 
La solución de dicha ecuación por el método del trapecio es: 
1 1 
=v +(—i + —1)h= 
Va+1 JA ( 2C La +1 IE al 


A AAN 


+1 
Vea 
Via E WEQ, FRE. Lao e 


Rea, , 1 = 2C 


Ecuaciones que dan como resultado el siguiente circuito equivalente in- 
cremental: 


Reqn+1 


n+1 Veq n 


Fig. 43. Circuito incremental del condensador. 
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Ya que los valores de corriente y tensión en el instante n son constantes 
para el cálculo de los valores en el instante n + 1. 

La resolución del problema requiere conocer los valores iniciales de 
tensión y corriente en el condensador para poder resolver el primer inter- 
valo. 


2. Bobina 
Igual que con el condensador hacemos con la bobina: 
A 
dt L 
+1 SE 1, id Gr +31 Ya h ER 
h h 
= G, 2L a 2L "+! 
h 
leq, =1 + =—v 
. n n 74 n 
a a E A AA h 


Fig. 44. Circuito incremental de la bobina. 


Del mismo modo que para el condensador, se necesitan las condicio- 
nes iniciales de estado de la bobina para resolver el problema. 


3. Resistencia 


El circuito incremental de una resistencia es ella misma, ya que en su 
ecuación no se hallan involucradas derivadas de ningún tipo. 

La sustitución de los componentes por sus circuitos equivalentes dará 
como resultado una red que sólo poseerá resistencias y generadores, de fá- 
cil análisis por el método nodal. Los coeficientes se determinarán en fun- 
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ción de la aproximación anterior. Con esto se consigue un sistema de ecua- 
ciones con coeficientes reales, cuya solución son los valores de las incóg- 
nitas del circuito en el siguiente instante de tiempo. Este proceso se repite 
cuantas veces sea necesario. 

Hay que resaltar la importancia que tiene el fijar un paso de integra- 
ción adecuado. Un paso demasiado grande nos llevará a errores en la so- 
lución, que harán que se aparte totalmente de la realidad. Un paso dema- 
siado pequeño, por otro lado, nos puede llevar a una duración del cálculo 
excesiva e innecesaria. 

Como ejemplo se puede ver en la figura 45 un circuito lineal y su cir- 
cuito incremental correspondiente, cuya obtención puede el lector justifi- 
car fácilmente. 


Fig. 45. Circuito lineal y su circuito incremental asociado. 


CIRCUITO LINEAL Y SU CIRCUITO 
INCREMENTAL ASOCIADO 


Los elementos no lineales se tratan del mismo modo que los lineales 
de cara al cálculo de la red incremental equivalente. La diferencia será 
que producirán una resistencia no lineal en ella, y, por tanto, el procedi- 
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miento de análisis tendrá que ser el visto en el capítulo «Análisis de cir- 
cuitos en continua», en la página 45. 

Tomemos como ejemplo el caso del diodo. Prescindiendo de efectos se- 
cundarios, lo trataremos como una resistencia no lineal, en la cual la re- 
lación corriente-tensión tendrá la siguiente forma: 


id, =1s. (e. /M- 1) 


n 


En la figura 46 se puede ver como ejemplo un circuito rectificador en 
media onda, con su circuito incremental equivalente. Dicho circuito se re- 
solverá en continua por el método iterativo, ya explicado con detalle en 
el apartado correspondiente. 


Rg + Vd - 


Fig. 46. Rectificador de media onda y circuito incremental asociado. 


En el caso de que haya transistores en la red, lo primero que se debe 
hacer es sustituirlos por los circuitos equivalentes que nos ofrecen los mo- 
delos vistos para el análisis en continua, formados por diodos y generado- 
res controlados, y proceder al análisis por el método general aquí estable- 
cido. 

En este caso es preciso considerar las capacidades presentadas en di- 
chos circuitos, que en el caso del análisis en continua se eliminaban. Como 
toda capacidad dará lugar a un circuito incremental equivalente como el 
de la figura 43. 

El problema es que estas capacidades no son lineales, por lo que la re- 
sistencia a la que dan lugar en el circuito incremental tampoco lo será. 
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Fig. 47. Organigrama del proceso de análisis transitorio de un circuito no lineal. 
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El tratamiento del problema es en estos casos bastante más complica- 
do. Normalmente, para una primera aproximación, los efectos capacitivos 
no se tienen en cuenta, usándose sólo cuando se quieren obtener unos re- 
sultados que se aproximen al máximo a la realidad. 

El lector interesado puede encontrar en el apéndice «Estudio de la ca- 
pacidad de una unión p-n», en la página 133, un estudio sobre la interpre- 
tación física de estas capacidades y el modo de tratarlas. 

En la figura 47 puede verse, como colofón de este capítulo, un organi- 
grama que resume el método de análisis transitorio de un circuito no li- 
neal. 
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A 


ASTA ahora se han tratado los métodos básicos para ana- 
lizar un circuito. El empleo del ordenador para realizar 
esta tarea resulta de gran utilidad, ya que permite abor- 
dar mayor complejidad y mayor volumen con el mismo es- 
fuerzo. 

Pero la solución a la que llegábamos no nos aportaba 
nada nuevo respecto a lo que suministraba el análisis tra- 
dicional. Gracias al ordenador podemos abarcar nuevos 
campos, como es el que se presenta en este capítulo, y que 
denominamos genéricamente sensibilidad. 

Buscamos con este tipo de análisis poder predecir el comportamiento 
del circuito ante cambios en uno o más componentes del mismo. Este es- 
tudio persigue una doble finalidad: 


q A -IIA4=>>=>AÁ 
o —K— 
q ___O_AA A AAA 
q$__E-dAázm AA AAA 
x$_x_____=z_——._€_[€[._.q_ nn —K—<— 
€ _ _ > -E-AAA<—]| 
_ÓÓÁÓÁ<Áá<AXfXf-z  Á€AKÁKÁA 
q$qEE____EAIH!HX REA XX 
q$E_A qm 
q _--_— És 
¡€  >_ ká e 
XK —— 
 _ ______-A«AAA<áá=>>. 
O. qq KA 


e Prevención de cambios indeseados. Un circuito, durante el transcur- 
so de su vida, puede ver variadas sus constantes por envejecimiento, efec- 
tos de radiación, cambios de temperatura, etc. Estas variaciones pueden 
hacer que se salga fuera de especificaciones, lo cual es una circunstancia 
indeseable en cualquier caso. 

Por otro lado, un circuito que se va a fabricar en serie debe estar dise- 
ñado considerando que no todos los componentes que se monten en cada 
unidad van a ser iguales. Hay unas tolerancias de fabricación contra las 
que debemos estar prevenidos, para evitar que un elevado tanto por cien- 
to de la producción resulte inútil. 

e Introducción de cambios deseados. Partiendo de un circuito base nos 
podemos preguntar cómo lograr unas mejores especificaciones en sus ca- 
racterísticas. Conociendo la sensibilidad de los parámetros a mejorar res- 
pecto a las diversas variables que componen el circuito, podremos llegar 
a deducir cuáles deben ser éstos para lograr el adecuado cambio en aqué- 
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llos. El estudio sistemático de estos métodos es lo que persiguen las téc- 
nicas de optimización. 


DEFINICIONES 


Consideremos como salida de un circuito el parámetro (tensión o 
corriente) cuyo estudio nos interesa. 

Dicha salida se puede expresar como una función de los valores de sus 
parámetros internos (resistencias, bobinas, condensadores...). Llamaremos 
a estos parámetros internos de un modo genérico con la letra «o», según 
lo cual: 


F(o1, 02, ... on) 

Su sensibilidad respecto a un determinado parámetro se define como 
el coeficiente que relaciona sus variaciones con la variación del mencio- 
nado parámetro. 

Llamando F a la variable salida, y ci al parámetro que varía, podemos 
escribir: 


AF = SF. Aoi 


Matemáticamente es, por tanto, la derivada parcial de la salida respec- 
to al parámetro en cuestión. 

A veces se suele utilizar en lugar de ésta la sensibilidad relativa, que re- 
laciona variaciones porcentuales en vez de variaciones absolutas. La re- 
presentaremos mediante una «s» minúscula. Se relaciona con la sensibili- 
dad absoluta a través de la expresión: 

e AF/E oi SE ai 
ei Aoi/oi F * óoi pom 

En el caso particular de que la función F sea una función de la red en 
el dominio de la frecuencia, podemos ponerla en forma módulo-argumen- 
to de la siguiente forma: 


Fo) = [FGo)| . e? 

En donde 0 será el argumento de F(¡w). La sensibilidad absoluta la po- 
dremos expresar como: 
Oz 
Sci 


Sci 


SF = 
E Sci 


SF, = el, SIFI +3 |F| elo, so 
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Y la sensibilidad relativa: 


=P +j0s> 
Llamaremos red adjunta o transpuesta de una red dada a aquella cuya 
matriz de admitancias es la transpuesta de dicha matriz en la red original, 
es decir, la misma en la que se han sustituido filas por columnas. 
Denominaremos a los parámetros de la red transpuesta con los mismos 
nombres que los de la red primitiva añadiéndoles «'». Según esto: 


ja Ty Y 
Daremos sin demostración la siguiente metodología de cálculo de las 
sensibilidades: 
Supongamos que tenemos un circuito con valor de salida Vo y quere- 


mos obtener su sensibilidad respecto a las variaciones de una admitancia 
Ya. 


1. Se analiza el circuito de interés, obteniendo el valor de tensión y 
corriente en bornas de la admitancia Ya, a los que llamaremos Va e la. 


2. Se analiza la red adjunta, con las siguientes modificaciones: 


e Sustitución de los generadores independientes de corriente por cir- 
cuitos abiertos. 

e Sustitución de los generadores independientes de tensión por corto- 
circuitos. 

e Introducción de un generador de corriente de valor unidad en el 
nodo de salida. 


Con esto, el sistema a resolver tendrá como matriz de admitancias la 
transpuesta del sistema original, y como vector de generadores equivalen- 
tes 
O, si el nudo i no es el de salida 


[leg], legi= 1, si el nudo i es el de salida. 


Nos interesa calcular en esta red también corriente y tensión en la ad- 
mitancia Ya. A estos valores les llamaremos V'a e Pa. 
3. La sensibilidad se calculará como: 
Vo pa q: , 
Sy, =-Va. Va 
En caso de que deseásemos calcular la sensibilidad de Vo respecto a 
una impedancia Za, el proceso sería el mismo, pero la fórmula final para 
la sensibilidad sería: 
Sy? = la. Pa 


En principio vemos cómo el problema del cálculo de sensibilidades se 
ha reducido a la resolución de dos redes. Esta doble resolución, «a prio- 
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ri», requeriría resolver dos sistemas de ecuaciones y, por tanto, el doble 
de trabajo de ordenador. 

La ventaja de este método se basa en que no es necesario realizar todo 
el proceso de solución del sistema de ecuaciones que genera la red adjun- 
ta. Si empleamos el método de descomposición L-U (ver «Resolución de 
sistemas de ecuaciones», en la página 133) tenemos: 


[Y]. [Y] = (1), 
que resolvemos como: 
[L] . [U] . [V] = [0] 


Ahora bien: recordemos que [Y”] es la matriz transpuesta de la [Y]. Por 
tanto: 


[Y] = (YT = ((L] . [UD” = [UT".. [LJ 
El sistema que tenemos que resolver será: 
[IV T=UL7, 
que se resolverá como: 
UT". [LJ".. [V'] = [1] 


El proceso se ha simplificado, ya que no es preciso descomponer la ma- 
triz de la red adjunta teniendo descompuesta la de la red principal. Será 
necesario realizar tan sólo la marcha inversa: 


[UT . [y] = [PJ 
[LJ”.. [V'] = [y] 


EJEMPLO DEL METODO DE CALCULO 
DE SENSIBILIDADES 


Para ilustrar el método de cálculo se va a resolver un sencillo ejemplo. 
Tenemos el circuito de la figura: 


Fig. 48. 


xi 


Deseamos conocer la sensibilidad de la tensión de salida Vo en función 
de las variaciones de la resistencia R1. Resolviendo directamente el sens 
cillo problema tendremos: 


R1 + R2 


R1 + R2 R1 + R2 
La sensibilidad la podemos calcular directamente a partir de su defini- 
ción en este sencillo caso: 


Vo = Vg 


SVo _ Vg R2 


RT RIO (RI+R2) 


Probemos ahora con la red adjunta o transpuesta. El circuito a resol- 
ver, según las reglas establecidas, será el siguiente: 


E Si 


Fig. 49. 


Nos interesa conocer de este circuito la corriente que atraviesa la re- 
sistencia Rl: 
1 ¿R2Z 
R1 + R2 


Y a continuación podemos obtener la sensibilidad respecto a R1 de la 
forma: 


Prl = 


0. A AS | 1 
RI+R2  RI+R2  (R1+R2) 


Que coincide con la que habíamos encontrado por el método directo. 

Basándose en los principios establecidos en este capítulo, la tarea de 
realizar un programa para cálculo de sensibilidades utilizando los progra- 
mas hasta ahora presentados no debe presentar ningún problema al lec- 
tor. Por tanto, se deja a su imaginación esta tarea, que puede servir como 
ejercicio para evaluar si realmente los conocimientos presentados se han 
adquirido de forma razonada y duradera. 


Sy =-Irl. rl =- 


OPTIMIZACION DE A 


NA vez que tenemos las herramientas básicas de análisis 
de cirucitos, vamos a abordar de forma general el proble- 
ma de la síntesis. El organigrama general de este proceso 
ya se presentó en la figura 2. En este capítulo nos vamos 
a centrar en el proceso iterativo que se debe llevar a cabo 
una vez definidos unos valores de componentes y analiza- 
do el circuito. Son varias las cuestiones que debemos res- 
ponder llegados a ese punto del diseño: 


A AAA e ¿Qué criterio se debe seguir para determinar si el 
circuito obtenido en el último análisis satisface nuestras exigencias o no? 
e ¿Qué componentes se deben alterar para acercarnos a las especifica- 
ciones buscadas? 
e ¿Cuál debe ser la magnitud y el sentido de ese cambio? 


En el diseño tradicional, estas preguntas se respondían en base a la ex- 
periencia del diseñador, que «intuía» el cambio necesario para hacer fun- 
cionar el equipo de acuerdo a lo establecido «a priori». La complejidad de 
cualquier método de cálculo hacía imposible abordar el problema de for- 
ma sistemática. 

El ordenador soluciona este problema, y en torno a él han surgido mé- 
todos para lograr realizar estos pasos de forma automática. El conjunto de 
estos métodos es lo que se conoce como técnicas de optimización de circui- 
tos. 

El diseño de acuerdo a estas técnicas tiene tres pasos fundamentales: 


1. Definir una función de comportamiento que sea una medida gene- 
ral del error entre la respuesta actual y la deseada de la red. 

2. Definir una red inicial, que será el punto de partida para el proce- 
so de diseño. 

3. Analizar la red y ajustar los parámetros de diseño deseados y la es- 
tructura de forma que se minimice la función de error. Este proceso de- 


ES 


berá repetirse hasta que la función de error se haya reducido por debajo 
de un cierto nivel o bien hasta que no se pueda reducir más. 


En este proceso vemos la aparición de dos nuevos problemas que aún 
no hemos tratado: 


e Cómo determinar la función de comportamiento o de error de la red. 
e Qué métodos emplear para encontrar sus valores mínimos. 


Estos dos temas son los que vamos a tratar en el resto del capítulo. 


DETERMINACION DE LA FUNCION DE 
COMPORTAMIENTO O DE ERROR 


La determinación apropiada de la función de error es fundamental en 
todo problema de optimización de circuitos. Si está mal definida, o no tie- 
ne en cuenta todas las variables posibles, el resto del procedimiento nos 
conducirá a unos resultados carentes de sentido físico. 

Con el fin de simplificar la notación a emplear llamaremos «o» a un vec- 
tor que contendrá los parámetros del diseño, es decir, los valores actuales 
de los componentes. Será un vector en columna, a efectos del cálculo ma- 
tricial, esto es: 


ol 
O= 02 sm (01, 02, 08) 
on 


La función de error más usada en los problemas de optimización es la 
de mínimos cuadrados Puede expresarse de forma general como: 


E(0) = (iy !W(0) . e(0,9)1? do (1) 


Expliquemos detalladamente cada miembro de esta igualdad: 


e (es una variable sobre la cual se optimiza el circuito. Puede ser la 
frecuencia, el tiempo, la temperatura... 

e e(o,p) es el error entre la respuesta actual y la deseada de la red. 
Como se ve, es función del vector de componentes y de la variable sobre 
la que actuamos. 

e W(() es una función de peso. Con ella queremos representar que no 
todos los valores de la variable q nos interesan de igual manera a la hora 
de calcular el error. Así, habrá unas zonas de mayor importancia, que lle- 
varán un valor de la función de peso mayor, y otras de importancia me- 
nor, en las que W(q) será de menor valor. 
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e [H) es el campo de variación de la variable q, es decir, el rango en 
el cual nos interesa optimizar la función. 


e La integral a lo largo de todo el campo nos indica que vamos a su- 
mar todos los errores en esa zona, y que ese total es lo que buscamos que 
sea mínimo. Se emplea una función cuadrática para evitar que al realizar 
esta integración los errores negativos se compensen con los positivos, dan- 
do como resultado un error nulo falso. 


Para aclarar mejor el sentido de cada uno de los términos de esta fun- 
ción de error vamos a estudiar un ejemplo de forma gráfica. Se trata de 
un filtro paso alto. 


Fig. 50. Función de error para un filtro paso alto. 


La figura (a) presenta las especificaciones de diseño. Se desea tener una 
atenuación de 100 en la banda eliminada y de 10 en la banda de paso. Es- 
tas especificaciones se deben cumplir hasta una frecuencia fmáx. 


US 


De lo dicho hasta ahora ya podemos decir, por pura analogía, cuáles 
son en este caso dos de los elementos componentes de la función de error: 


e La variable q va a ser la frecuencia, que representaremos como f. 
e El rango (H) es el intervalo de frecuencias de O a fmáx. 


En la figura (b) se ha superpuesto a la función que representaba las es- 
pecificaciones de diseño la función obtenida en la primera aproximación. 
Como vemos, no coinciden, y en la figura (c) se representa la diferencia 
entre ellas, que constituye el error, designado en (1) por e(o,p). 

Continuemos analizando el ejemplo. El error que más nos interesa mi- 
nimizar en este caso es el cometido en la banda pasante, ya que la señal 
que se vea afectada por él es señal útil de la salida del circuito. Introduci- 
mos para ello una función de peso, W(f), representada en la figura (d), me- 
diante la cual triplicamos el error para f > fc. En (e) se ve el efecto mul- 
tiplicador señalado. 

Por último, al elevar al cuadrado convertimos todos los errores en po- 
sitivos, y obtenemos como resultado de la integral el total del área rayada 
en (f). 

Normalmente el rango de valores permitidos de la variable indepen- 
diente es discreto en lugar de continuo. En este caso la integral se con- 
vierte en una suma, y tendremos: 

E(o)= E |W,.e(o)l? (1) 
ic (HI) 

A esta función se le denomina función objetivo, ya que su mínimo nos 
dará la mejor solución posible. ¿ representa cada uno de los puntos de la 
variable y que debemos considerar para la optimización. 


METODOS PARA HALLAR EL MINIMO 
DE LA FUNCION OBJETIVA 


Casi todos los métodos numéricos para determinar el óptimo de una de- 
terminada función son iterativos. Parten de una solución inicial dada y, a 
partir de ella, van generando una serie de estimaciones de dicho vector. 
Cada una de ellas debe presentar una mejora respecto a la anterior en lo 
que se refiere al valor de la función objetivo que proporciona. 

Los distintos métodos difieren entre sí en la estrategia para producir la 
serie de estimaciones. Podemos establecer dos tipos fundamentales de mé- 
todos: 


+ Métodos de búsqueda directa: Se basan únicamente en comparación 
de valores de la función objetivo. Suelen ser simples y fáciles de progra- 
mar, pero lentos. 


e Métodos de gradiente: Usan las derivadas de la función objetivo para 
lograr una convergencia más rápida. Los hay de primero y segundo orden, 
según que utilicen sólo las primeras derivadas o las primeras y las segun- 
das, respectivamente. 


METODOS DE BUSQUEDA DIRECTA 


En ellos se trata de reducir el valor de la función objetivo E usando 
comparaciones del valor que toma en una estimación con los valores en 
puntos próximos. Esto es: 


e Dado un componente del circuito, se estudia en qué direccion debe 
variar su valor (aumentar o disminuir) para producir un descenso en la 
función objetivo. 

e Este proceso se repite para cada uno de los componentes, obtenien- 
do así una «dirección global de búsqueda». 

e Una vez obtenida esta dirección se realiza un paso en ella, esto es, se 
actualizan los valores de los componentes en la dirección elegida y se cal- 
cula la nueva función objetivo. 

e Estos pasos se repiten hasta que las variaciones de la función objeti- 
vo que se logren en cada paso sean menores que un valor máximo de error 
fijado con antelación. 


Hay dos formas básicas de plasmar este procedimiento: 
e Mediante una minimización independiente para cada parámetro. 


e Se comienza con una variable y se prosigue con ella hasta que sus 
variaciones en ambos sentidos aumenten la función de error. 

e Este proceso se repite para el resto de las variables. 

e Una vez realizado con todas se vuelve a repetir el ciclo completo has- 
ta que al terminar uno de ellos no se haya conseguido una mejora 
por encima del margen de error. 


Este ajuste puede ser, por tanto, muy complejo si existen interacciones 
entre los distintos elementos de cara a la función objetivo. 


e Por el procedimiento de variables alternadas. En lugar de optimizar 
al máximo para cada variable, se realiza un paso en una de ellas y se pasa 
a evaluar la dirección a continuar en la siguiente. Este proceso se repite 
ciclicamente hasta llegar a un valor de la función de error imposible de 
minimizar. 


Hay otros muchos métodos de optimización por búsqueda directa, para 
cuyo estudio se remite al lector a la bibliografía especializada. 


1:1:7 


METODOS DE GRADIENTE 


Si bien estos métodos conducen a unos resultados aceptables de una 
forma más rápida, son mucho más difíciles en su concepción, por lo que 
aquí se van a tratar de forma muy «elemental. 

Estos métodos se basan en el desarrollo de la función de error del cir- 
cuito, que hemos denominado e(o), en serie de Taylor en un entorno de 
la solución actual. Con eso podremos obtener: 


e(o + 40) = elo) +8T.A0+—=- Ac".[H].o0 


en donde: 
de de Se 
T = 
a E 
Se Se Se 
801? 801802 eN $0180n 
Se Se Se 
(H] = 802801 802? dé 30280n 
Se Se Se 
¿ongo1  Sdondods2 Son? 


g es el gradiente o jacobiano de la función y [H] el Hessiano. Con este de- 
sarrollo tenemos: 


e Métodos de primer orden: Toman sólo los dos primeros términos, 
sin considerar el que contiene la matriz [H]. El más conocido es el llama- 
do método del descenso más rápido. 

En este método, a partir de una solución dada oi, se obtiene la siguien- 
te restando a ésta el producto de una constante por el valor del vector gra- 
diente en el punto oi: 


0; ,1=0,-K;. g(0,) 


La constante puede variar de una iteración a otra, y se determinará de 
forma que o + 1 dé como resultado una función objetivo menor que oi. 

El cálculo de las derivadas que componen la matriz jacobiana lo pode- 
mos realizar de muy diferentes maneras. Para ilustrar la más cómoda para 
nosotros y que nos permite aplicar conocimientos que ya hemos adquiri- 
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do, reconsideremos el significado que la función e(o) tiene en nuestro cir- 
cuito. Es la diferencia entre nuestra salida y la salida que deseamos, esto es: 


e(o) = V (0) - Vop 


En donde V puede ser cualquier función de la red (tensión o corriente) en 
un nodo o en un determinado componente. Vop es constante respecto a 
los parámetros de diseño. Por tanto, al calcular la derivada parcial tendre- 
mos: 

Se 


la) 
e Va) > - 
S0i Soi yde) = epa 


SV(o) 


= Sy, 
Sci o 


Expresión de la sensibilidad de la función de salida respecto al parámetro 
de la red oi. 

Para calcular el jacobiano de la función bastará por tanto hallar las sen- 
sibilidades de la red respecto a todos los componentes que la forman, que 
se pueden obtener fácilmente por el método de la red adjunta. 


e Métodos de segundo orden: toman el desarrollo en serie de Taylor in- 
cluyendo el Hessiano. De los múltiples algoritmos existentes vamos a ex- 
poner la metodología de cálculo de uno de los más conocidos y más efi- 
cientes: el de Davidon-Fletcher-Powell. 

Este método no utiliza puramente la matriz Hessiana, sino que emplea 
una matriz que llamaremos ¿Sñ, que se actualiza en cada iteración. El pro- 
ceso para la iteración [1] es el siguiente: 


1. Se calcula la dirección de búsqueda haciendo 
d, =-1S]. g(0) 


En la primera iteración se toma para [S] la matriz identidad, con lo que 
la dirección de búsqueda coincide con la de los métodos de primer grado. 

2. Se efectúa una búsqueda del mínimo de la función objetivo en esa 
dirección. Esto se dará para un valor: 


5,10 PK .d,=0,+p, 

En otras palabras, buscamos un valor de K tal que E(ci + 1) sea el mí- 
nimo en la dirección di. Este mínimo se calcula por búsqueda directa, dan- 
do valores a K y calculando el valor de la función objetivo. 

3. Calcular el nuevo jacobiano de la función en este nuevo punto, 
g(oi + 1). 

4. Obtener la diferencia entre los jacobianos de las dos iteraciones: 


9; e glo, 1) a g(o;) 
5. Calcular la nueva matriz ¿Sñ mediante la fórmula: 
[S;, ,] = [S¡] + [4] + [B,), 
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en donde 


P; - Pi _ [S].q q -[S] 
Y A AAA 
Pp; -9 q; [Si]. 9, 
La complejidad de estas fórmulas está más en la forma de escribirlas 
que en el cálculo necesario para obtenerlas. Notar, por ejemplo, que los 
denominadores que aparecen en las expresiones de las matrices [Ai] y [Bi] 


son números, no matrices, con lo cual no hay que realizar ninguna inver- 
sión. Por ejemplo, en la matriz [Ai]: 


P-9,=P,-9,+P,-9,+ +. +P,- 9, 


6. Este procedimiento se repetirá desde el paso (1) hasta que cada 
componente de di sea menor que una cantidad prefijada. 


[A] = 


El hecho de que este método sea de segundo orden es porque la serie 
de matrices [Sk] que calculamos tiene por límite el Hessiano de la fun- 
ción. La interpretación de dicha serie es la siguiente: 


e Para las primeras iteraciones, ¿Sñ es la matriz identidad, con lo cual 
la búsqueda se hace por el método del descenso más rápido. Este método 
es bueno cuando se está relativamente lejos de la solución, pero es lento 
en proximidades de dicho punto. Por eso: 

e En las iteraciones finales, [S] se aproxima al valor del Hessiano de: 
la función. Estamos empleando un método de segundo orden, muy bueno 
en las proximidades de la solución, pero con problemas para valores leja- 
nos. Este sistema aprovecha lo mejor de los métodos de primer y segundo 
orden. 


CONCLUSIONES 


Este capítulo no es sino un brevísimo repaso de las técnicas de optimi- 
zación de circuitos para dar una idea de su complejidad y de los proble- 
mas que plantean. Es cierto que están más cerca que ningún otro sistema 
de la automatización total del diseño, pero necesitan una interacción muy 
fuerte con el usuario, que es quien debe determinar: 


e Estructura de la red y primera solución. 
e La función objetivo más adecuada para resolver el problema. 
e El algoritmo de búsqueda del mínimo utilizado. 


Una vez establecido esto no ha terminado el trabajo, pues muchas ve- 
ces los métodos empleados presentan difíciles problemas de convergen- 
cia, estabilidad, etc., que pueden hacer necesaria su sustitución a partir de 
una determinada iteración. 
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APENDICE A. ESTUDIOS 
ESPECIALES == 


ANALISIS NODAL DE UN CIRCUITO 


Se considera, para simplificar el desarrollo, que el circuito a analizar 
esté compuesto de los siguientes elementos: 


Resistencias. 

Bobinas. 

Condensadores. 

Fuentes de corriente independientes. 

Fuentes de corriente controladas por tensión. 


El circuito consta de m ramas y n nudos. 
Una rama generalizada tendrá la siguiente estructura: 


Fig. 51. Rama generalizada. Convenio de signos. 


1:21 


La rama generalizada consta, por tanto, de 


+ Un generador independiente de corriente. 
e Un generador independiente de tensión. 
+ Una admitancia pasiva (o un generador de corriente dependiente). 


Las ecuaciones básicas de esta rama, considerada como la número h» 
del circuito, serán: 
eh = Vh - Eh 
ih = —Ih + jh = -Ih + Yh (eb + Eb) + gm (el + El) (1) 


A continuación se definen las siguientes matrices: 


[Igen]: Generadores de corriente. Dimensión: m X 1. Cada Igeni represen- 
ta el valor del generador de corriente asociado a la rama i. 
[Egen]: Generadores de tensión. Igual dimensión y significado que ¿Igenñ. 


[ele Vector de tensiones en cada rama. Dimensión: m xX 1. 
[1]: Vector de corrientes en cada rama. Dimensión: m x 1. 
[Y] Matriz de admitancias de rama. Dimensión: m x m. Un elemento 


de esta matriz Yijtendrá por valor 


e Yi, sii=j (Admitancia de la rama i). 
e 0)sii<>j. 
[gm]: Matriz de conductancias de los generadores de corriente contro- 


lados por tensiones de rama. Dimensión: m x m. Un elemento ge- 
nérico gmij tendrá por valor: 


e gmi, si el generador de corriente se halla en la rama i y está con- 
trolado por la tensión en la rama j. 
e (, en caso contrario. 


[Vn]: Vector de tensiones de nudo. Dimensión n x 1. 
[A]: Matriz de incidencia. Representa la geometría del circuito. Dimen- 
sión: n x m. Sus elementos Aij tienen como valor: 
e 0, si la rama j no está conectada al nudo i. 
e 1, si la corriente de la rama j sale del nudo i. 
e -1, si la corriente de la rama j entra al nudo i. 


De las definiciones de las matrices se pueden establecer las siguientes 
relaciones: 


[e] = [AT'.. [Vn] (2) 
[A]. [1] =0, (3) 


según la ley de Kirchoff de corrientes. 
Por otra parte, el sistema de ecuaciones de corrientes de (1) de todas 
las ramas del circuito se puede expresar en forma matricial como: 


[1] = -[Igen] + [Y] . [e] + [Y]. [Egen] + [gm] . [e] + [gm]. [Egen] — (4) 
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considerando la existencia de admitancias y generadores controlados si- 
multáneamente. Sustituyendo (4) en (3) se obtiene: 


-[A] . [Igen] + [A] . [Y] . [e] + [A] . [Y] . [Egen] + [A] . [gm] . [e] + 
+ [A]. [gm] . [Egen] = 0 


Transformando esta ecuación: 
[A] . [Y+gm]. [e] = -[A] . [Y+gm] . [Egen] + [A] . [Igen] 
Sustituyendo (2) en esta última ecuación obtenemos: 
[A]. [Y+gm]. [AJ" . [Vn] = [A] . [Y+gm] . [Egen] + [A] . [Igen] 
Ecuación que se puede presentar de la forma: 


[Yn] . [Vn] = [eq] (5) 
En donde 
[Yn] = [A] . [Y + gm]. [AJ" 


es la matriz de admitancia nodal reducida (dimensión n x n) 
[leg] = -[A] . [Y+gm]. [Egen] + [A] . [Igen] 
es el vector de generadores de corriente equivalentes (dimensión n x 1) 


Las ecuaciones representadas en notación matricial en (5) son las ecua- 
ciones nodales de un circuito. 

Para aclarar el procedimiento a seguir se va a realizar un ejemplo. El 
circuito de la figura es el equivalente en pequeña señal de un transistor co- 
nectado en emisor común: 


Fig. 52. Ejemplo de cálculo. 
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El circuito consta de cinco ramas y tres nudos. En la figura, las ramas 
se han indicado con numeración romana. Las matrices de interés son las 
siguientes: 


e Matriz de generadores equivalentes: Sólo posee un término, ya que 
sólo hay un generador independiente: 


lg 
[leg] = 0 
e Matriz de incidencia: 
OO 
[A] = 0 -1 1 -1 
0-0 A dl 
e Matriz de admitancias: 
1/Rg 0 0 0 0 
O 1/hie O 0 0 
[Y] = 0 O 1/Re O 0 
0 0 0 0 0 
0 0 0 O. 1/Re 
e Matriz de transconductancias: 
0 0.0 0 
0 0 0" 0 
[gm] = 0 O 0 0 
0 gm 0 0 0 


0.0 0.0.0 
Cálculo de la matriz de admitancia nodal reducida: 
[Yn] = [A]. [Y + gm]. [AJ' = 


l/Rg 0 0 0 0 1 00 

1.10 00 0 1/hie 0 0 0 10 
MAA 0 0 1/Re 0 0 010 = 

0.00 11 0 gm 0 0 0 Meri 

00 0. 0 1k 0 01 
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1 0 0 
1/Rg 1/hie 0 0 0 1 -1.0 
= O -1/hie-gm 1/Re 0 0 0 1.0 
0 gm 0 O 1/Re 0 -1 1 
0 o. 1 
1/Rg + 1/hie -1/hie 0 
= -1/hie - gm 1/hie+ 1/Re+gm 0 
gm -gm 1/Re 


Esta matriz se puede descomponer en otras cuatro que muestran la 
aportación de cada una de las ramas: 


1/Rg 0 O 1/hie -1/hie 0 
[Yn]= 0 0 0 + -1/hie 1/hie 0 + 
0 0 0 0 0 0 
0 0 0 0.0 0 
+ -gm gm 0 + 0: 0-0 
gm -gm 0 0 0 1/Re 


Que responden a las reglas establecidas en el capítulo «Fundamentos del 
diseño de circuitos asistido por ordenador», en la página 13. 


El sistema de ecuaciones obtenido es, por tanto: 


1/Rg + 1/hie -1/hie 0 VI lg 
-1/hie - gm 1/hie + 1/Re+gm 0 v2 = 
gm -gm 1/Re V3 


O bien en la forma habitual: 


(1/Rg + 1/hie) V1 - (1/hie) V2 = lg 
(-1/hie - gm) V1 + (1/hie + 1/Re + gm) V2=0 
gm V1 - gm V2 + (1/Rc) V3 =0 


RESOLUCION DE UNA ECUACION DIFERENCIAL 
POR EL METODO DEL TRAPECIO 
Se desea resolver la ecuación diferencial en forma explícita: 


y' = f(x) 


EZ 


De la cual conocemos un punto, que llamaremos (X0, YO). Este punto 
es la condición inicial del problema. 

Buscamos obtener el valor de la función y para un determinado con- 
junto finito de valores de x (x1, x2, x3...) equidistantes entre sí, de tal ma- 
nera que: 


x2-x1=x3-x2=x4-x3=..=h 


h recibe el nombre de paso de integración. 

Como sabemos, la derivada de una función nos representa la pendien- 
te de la recta tangente a ella en un determinado punto. Nuestra ecuación 
diferencial nos da, por tanto, el valor de las tangentes a la función que bus- 
camos en cada uno de sus puntos. 

Todos los métodos de resolución numérica de ecuaciones diferencia- 
les se basan en el mismo principio: sustituir la función que no conocemos 
por una recta cuya pendiente sea el promedio de las pendientes de las tan- 
gentes a la curva desconocida en distintos puntos internos al intervalo 
(xi, xi + h). 

En el caso del método del trapecio, este promedio es simple de calcu- 
lar, y es igual a la media aritmética de los valores de las derivadas en los 
puntos extremos del intervalo. Matemáticamente: 


y )+y 0,1) 
2 


La interpretación geométrica de este método la podemos ver en la fi- 
gura 53. 


O (X;, 1 X) 


Y 


Y (xi) 


Xx; X,+h 


Fig. 53. Interpretación gráfica del método del trapecio. 
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El punto yi + 1 es la solución que nos da nuestro método. Si la curva 
real fuese la dibujada, tendríamos un error en esta aproximación de valor 
Ei + 1, que además se iría acumulando en el siguiente intervalo. 


Como hemos fijado que el paso de integración sea constante, podemos 
poner la expresión del método de la forma: 


Vis: y, + 10 A h 


ESTUDIO DE LA CAPACIDAD DE UNA UNION P-N 


En todo dispositivo de unión p-n existen unas acumulaciones de por- 
tadores en su volumen, irregularmente distribuidos. Estas acumulaciones 
varían de forma no lineal con la tensión que se aplique al dispositivo. 

El fenómeno se manifiesta al exterior como un efecto capacitivo. Pues- 
to que la variación de carga no es lineal, la capacidad equivalente que pre- 
senta el dispositivo tampoco lo será. 

Llamemos Cpn a esta capacidad. Se puede considerar constituida por 
otras dos básicas, fruto de dos fenómenos físicos distintos: 


Ct = Cd + Cz 


e Cd: Capacidad de difusión. Se debe a la existencia de portadores de 
corriente en el volumen del componente, lejos de la zona de la unión. Este 
es el efecto fundamental bajo polarización directa. La expresión para ella 
s: 


9) 


Tt 
Cd = —1s (ev 1 
ye (e ) 


donde tt es el tiempo de tránsito de los portadores. 


e Cz: Capacidad en la zona de carga espacial o zona de la unión pn. Su 
importancia relativa es mayor bajo polarización inversa. Se relaciona con 
la tensión en los extremos del diodo a través de la expresión: 


K 


C= Ya, -v 


Q., es la llamada tensión de banda plana del dispositivo y K una constante 
que depende de la geometría del componente. 
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El circuito incremental de una capacidad no lineal es exactamente el 
mismo que el del caso lineal, excepto que la resistencia de valor constante 
obtenida allí será no lineal: 

=V, + - i 
2) + 


h 
A 
¿+1 2C(Vn + 1) 


Veg 


Para analizar este componente deberemos sustituirlo por una rama li- 
nealizada, según se estudió en el apartado «Análisis de circuitos en conti- 
nua», en la página 45. 


RESOLUCION DE SISTEMAS DE ECUACIONES 


Muchos son los métodos para resolver sistemas de ecuaciones. En este 
capítulo sólo se va a tratar uno de ellos, debido a su uso generalizado en 
los programas de análisis de circuitos. Este método es el llamado de Cho- 
leski o de la descomposición L-U. 

Su fundamento es muy sencillo. Tenemos el sistema: 


[A]. [X] = [B] 


Que es de orden n. Supongamos que podemos descomponer la matriz cua- 
drada [A], de orden n x n en el producto de otras dos matrices: 


[A] = [LJ]. (U] 


Si esto es posible, el sistema de ecuaciones que tenemos lo podemos 
representar por: 


(L] . [U] . (X] = [B] 


Basándonos en esto, podemos descomponer el sistema de ecuaciones 
en otros dos, de forma que: 


[0]. [(X] = [Y] (1) 
[E]. [Y] = [B] (2) 


Con esto, hemos descompuesto el problema en la resolución de dos siste- 
mas. Primero habrá que resolver el (2) y, con el vector [T] que allí se ob- 
tiene, se resolverá el sistema (1). 

¿Dónde está la ventaja? Hasta ahora lo único que se ha conseguido es 
duplicar aparentemente el problema. 
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La gran ventaja de esta descomposición es que se buscan las matrices 
[1] y [u] de forma que sean triangulares, esto es, de la forma: 


11 0 0 dd 0 
Pr 12 30 0 0 
[L] = 81 182 "B3 (0) 
Ini In2 1In3 .. Inn 
i u2 u3 . uln 
0 123. mula 
== 0 0 1 Ba 
0) 0 0 se 1 


Los sistemas de ecuaciones con matrices triangulares son de resolu- 
ción sencilla. Basta con despejar la incógnita de la ecuación en la que se 
encuentra sola para calcular su valor. Después, este valor calculado nos sir- 
ve para sustituir en la siguiente ecuación y así obtener una nueva incóg- 
nita. Este proceso se repite hasta resolver el sistema. 

Los coeficientes de las matrices [L] y [U] se obtienen directamente por 
igualación del producto de ambas matrices a la matriz [A]. Los resultados 
que se obtienen se pueden resumir en las siguientes fórmulas: 


1 = a, - A th» ji 
k=1 
1 i-1 
Y 07 2 by.» 1<j 


Este método se puede ver con mayor claridad si realizamos un sencillo 
programa para resolver sistemas de ecuaciones que lo utilice. La rutina de 
resolución que contiene, por otra parte, es la que se empleó en el progra- 
ma de resolución de circuitos en continua. 


E 


rr —— A A PP 
(* pste programa resuelve un sistema ae ecuaci: 
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fama = arravil1l GOGrdan Maxima 1 (Ordan Mawime,] 
cema — array. .Vvraen Max1moO0,1..vraen mMaXx1mo ] 


A 
elve S5istemaívar A: Sistema 


-A[i,k]*A[k,3]; 
end; 
for i:=j+1 to Orden do 
begin 
(* Coeficientes matriz U *) 
for k:=1 to j-1 do A[j,i]:=A[j,1i] 
AG KA 
A AAA 
end; 
end; 


(* Marcha inversa. Solución de LY = b *) 


For i:=1 to Orden do 
begin 
Y[i]:=B[i]; 
for k:=1 to i-1 do Y[ilJ:=Y[i] - Ali,k]*Y[k]; 
A NA Z 
end; 


(* Marcha inversa. Solución de UX = Y *) 


For i:=0rden downto 1 do 
begin 
X[lil:=Y[i]:; 
for k:=0rden downto i+1 do X[i]:=X[i] 
"AE XP 


(*====FIN RESOLUCION DE-UN SISTEMA DE ECUACIONES====*) 
* 


(* Variables globales del programa *) 


Var 
A: Sistema; (* Matriz del sistema *) 
B, (* Términos independientes *) 
X: Vector; (* Resultados *) 
N: integer; (*-—Orden del sistema *) 
i,j: integer; (* Variables auxiliares *) 
BEGIN 
WriteLn('RESOLUCION DE SISTEMAS DE ECUACIONES* y 
Writeln(!-=-==-====-== oo o-ooooooooo-o-- 3 
WriteLn; 
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Write('Orden del sistema...... >'); 
ReadLEn(N); 


WriteLn; 
WriteLn('Introduzca los coeficientes del sistema:'); 


For i:=1 to N do 
For j:=1 to N do 


begin 
LA A A AAA >); 
ReadLn(A[i,j]); 
end; 
WriteLn; 


NriteLn('Introduzca los términos independientes:'); 


For i:=1 to N do 


begin 
A A A AE == 
ReadLn(B[i]); 

end; 


Resuelve Sistema(A,B,N,X); 


WriteLn; 
WriteLn('Descomposición del sistema:'); 


WriteLn; 
WriteLn('Matriz L:'); 


for i:=1 to N do 
LA OA 
begin 
E e A PS | 
end; 


5 A11,31) 


WriteLn; 
WriteLn('Matriz U:'); 


for i:=1 to N do 
For j:=i+1 to N do 
begin 
IA A | 
end; 


*,A11,31); 


WriteLn; 
WriteLn('Solución del sistema:'); 
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Consideremos ahora el caso, muy frecuente en análisis de circuitos, de 
que tengamos un sistema de ecuaciones complejas, esto es, en el que las 
variables y sus coeficientes son números complejos. Tenemos dos proce- 
dimientos para atacarlo: 


e Descomponer el sistema en otro de doble orden, igualando en cada 
ecuación por separado parte real y parte imaginaria de ambos miembros. 

e Resolver el sistema de igual manera que el real, pero operando con 
números complejos. 


Esta última solución es la que se ha llevado a cabo en el programa que 
a continuación vamos a analizar, y cuya rutina principal se ha empleado 
al analizar los métodos de resolución de circuitos en alterna. 

Puesto que en Pascal estándar no existen las variables de tipo comple- 
jo, se ha definido un tipo de datos para sustituirlas, que consiste en una ma- 
triz de dos números reales. Al mismo tiempo se han realizado los proce- 
dimientos para ejecutar las cuatro operaciones básicas entre números 
complejos y el paso de rectangulares a polares. De este modo, la misma 
estructura del programa de resolución para sistemas reales sirve aquí, sin 
más que sustituir: 


e Los elementos de las matrices de reales, por otros del nuevo tipo 
efinido. 

e Las operaciones entre dichos elementos: De las normales en Pascal 
a las implementadas en los nuevos procedimientos programados. 


[em 


13 


W 


Orden Maximo = 20; 
(* Orden máximo permitido al sistema*) 


Type 
Complejo= array[1..2] of real; 
(* Variable compleja *) 
Sistema = array[1..Orden Maximo,1..Orden Maximo] 
of Complejo; 
(* Tipo de datos para la Matriz del sistema *) 


Vector = arrayll1..Orden Maximo] of complejo; 
(* Tipo de datos para los vectores de incógnitas, 
terminos independientes... *) 


Procedure R P(Var C:Complejo); 
(* Rectangulares a polares *) 
var 
R: real; 
begin 
R:=Sgrt(Sqr(C[1])+Sqr(C[2])); 
if Abs(Cc[1]) < 1E-10 


then 
begin 
AA 
then C[2]:=Pi/2 
else C[2]:=-Pi/2; 
end 
else 
begin 
if cf1]>0 
then C[2]:=ArcTan(C[2]/C[1]) 
else C[2]:=ArcTan(C[2]/C[1]) + Pi; 
end; 
Ct]: =R; 
end; 


Procedure IguC(C1: Complejo;var Cr:Complejo); 
(*—Iguala complejos *) 


var 
i:integer; 
Begin 
For i:=1 to 2 do Crli]:=C1[i]; 
end; 


Procedure SumCc(C1,C2: Complejo;var Cr:Complejo); 
(* Suma de complejos *) 
var 
i:integer; 
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Begin 
For-i:=to-2-—do Crtid:=CHki4c2t3F; 
end; 


Procedure ResC(C1,C2: Complejo;var Cr:Complejo); 
(* Resta de complejos *) 


var 

i:integer; 
Begin 

For i:=1 to 2 do Cr[i]:=C1[i]-C2[i]; 
end; 


Procedure MulC(C1,C2: Complejo;var Cr:Complejo); 
(* Producto de complejos *) 
var 
i:integer; 
Begin 
Cr[1]:=C11[1]*C2[1]-C1[2]*C2[2]; 
Crf2]:=C1[1]*C2[2]+C1[2]*C211]; 
end; 


Procedure DivC(C1,C2: Complejo;var Cr:Complejo); 
(* División de complejos *) 


var 
i:integer; 
D: real; 

Begin 
D:=sqr(C2[1])+Sqr(Cc2[2]); 
Ccr[1]:=(c1[1]*c2[1]+C1[2]*C2[2])/D; 
Cr[2]:=(-Cc1[1]*C2[2]+C1[2]*C2[1])/D; 


end; 


(* Procedimiento para resolver un sistema de 
ecuaciones mediante factorización L - U 
Parámetros: 
-=A : Matriz con los coeficientes del sistema. 
En ella se devuelve el 
resultado de la factorización. 


-B : Vector de términos independientes. 
-= Orden: Orden del sistema. 
-X : Vector de soluciones *) 


Procedure Resuelve Sistema(var A: Sistema; 
var B: vector; Orden: Integer;var X:Vector); 


var 


i,j,k : integer; 
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y Vector;z 
(* Solución intermedia de LY = B 
c : Complejo; 
Begin 


(* Marcha directa. Descomposición A = LU 


For j:=1 to Orden do 


Begin 
for i:=j to Orden do 
begin (* Coeficientes matriz L 
for k:=1 to j-1 do 
begin 
MulC(A[li,k],A[k,j]1,C); 
ResC(Ali,j],C,A[li,j3]); 
end; 
end; 
for i:=j+1 to Orden do 
begin (* Coeficientes matriz U 
for k:=1 to j-1 do 
begin 


MulC(A[j,k],A[k,i],C); 
ResC(A[j,i],C,A[j,i]); 


end; 
Divc(A[3,1],A13,3],A(3,1]); 
end; 
end; 


(* Marcha inversa. Solución de LY = b *) 


For i:=1 to Orden do 


begin 
IguC(B[i],Y[i]); 
for k:=1 to i-1 do 
begin 
MulC(A[i,k],Y[k],C):; 
ResctiEEXED; 
end; 
Divc(Y[i],A[i,i],Y[i]); 
end; 


(* Marcha inversa. Solución de UX = Y *) 


For i:=0rden downto 1 do 
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begin 
IguC(Y[i],X[i]); 
for k:=0rden downto i+1 do 
begin 


= 


5 


+) 


5 


MulC(A[i,k],X[k],C); 
ResC(X[i],C,X[i]); 
end; 
end; 
end; 


=============o======== o o=o===============_o=o======= 


(*====FIN RESOLUCION DE UN SISTEMA DE ECUACIONES====*) 
* 


(* Variables globales del programa *) 


Var 
A: Sistema; (* Matriz del sistema *) 
B, (* Términos independientes *) 
X: Vector; (* Resultados *) 
N: integer; (* Orden del sistema *) 
i,j: integer; (* Variables auxiliares *) 


Procedure ReadLnC(var C: Complejo); 
(* Lectura complejo 


begin 
Write('Parte Real: '); 
Read(C[1]); 
Write(' Parte Imaginaria: '); 
ReadEIn(C[2]); 

end; 


BEGIN 


WriteLn('SISTEMAS DE ECUACIONES COMPLEJAS*' = 
Wri teLn(' A 3 i 
WriteLn; 


Write('Orden del sistema...... a 
ReadLn(N); 


WriteLn; 
WriteLn('Introduzca los coeficientes del sistema: 


For i:=1 to N do 
For j:=1 to N do 
begin 
ARAU E WS >; 
ReadEInc(A[i,j]); 
end; 


137 


WriteLn; 
WriteLn('Introduzca los términos independientes:'); 


For i:=1 to N do 
begin 
NES PASA RS 3 
ReadEInC(B[i]); 
end; 


Resuelve Sistema(A,B,N,X):; 


WriteLn; 
WriteLn('Descomposición del sistema:'); 


WriteLn; 
WriteLn('Matriz L:'); 


for—3i:=1-to-N-do 
For j:=1 to ido 
begin 
A O A SR + 
11 


end; 


WriteLn; 
WriteLn('Matriz U:'); 


for-1i:=1-to N- do 
For j:=i+1 to N do 
begin 
Hriteln('0[',i,*',*,3,"] 2 


end; 


Writebn; 
WNriteLn('Solución del sistema:'); 


For i:=1 to N do 
A E A 0 SS OS 


end. 
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NOTAS PARA ADAPTAR LOS PROGRAMAS DEL 
PRESENTE LIBRO A SU ORDENADOR === 


Como ya se expuso al comienzo del libro, todos los programas se han 
realizado en lenguaje Pascal. Han sido compilados mediante TURBO PAS- 
CAL en un ordenador IBM-PC. 

Si usted dispone de un compilador de Pascal para su ordenador no en- 
contrará prácticamente ningún problema en ejecutar dichos programas, 
ya que incorporan instrucciones que son comunes a todos los compilado- 
res de Pascal existentes en el mercado. Unicamente puede necesitar agru- 
par las declaraciones de constantes, tipos y variables al comienzo del pro- 
grama. Para ello: 


e Busque dentro del programa todas las definiciones de constantes. Las 
reconocerá por ir precedidas por el identificador CONST y terminar con 
otra palabra reservada, que podrá ser BEGIN, TYPE, VAR o PROCEDURE. 
Tome, por tanto, todas estas líneas y colóquelas en el programa tras el iden- 
tificador de programa que se encuentra en la primera línea. 

e Haga lo mismo con las definiciones de tipos. Estas las encontrará tras 
la palabra reservada TYPE. Coloque dichas definiciones tras el bloque de 
definición de constantes antes agrupado. 

+ Repita este mismo procedimiento con las definiciones de variables 
globales, esto es, aquellas que no se encuentran introducidas dentro de 
un procedimiento. 


Con esta sencilla transformación podrá usted ejecutar su programa tras 
compilarlo con los compiladores de Pascal más comunes del mercado. 

Si, por el contrario, usted sólo posee un intérprete de BASIC, la tarea 
puede resultar un poco más compleja, pero en ningún caso imposible. 
A continuación se ofrece una serie de normas maestras para ejecutar dicha 
conversión: 


e Cambie los nombres de las variables al máximo número de letras que 
permita su ordenador. Este número es normalmente 2. 
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e Introduzca las variables alfabéticas con el identificador adecuado al 
intérprete BASIC que utilice. Normalmente deberá terminarlas con un «$». 
e Sentencias de asignación: Mientras en PASCAL se utiliza para asig- 
nar a una variable el símbolo «:=», en BASIC se utiliza solamente «=»: 


o —_  ___A<______———  —_—__— — _ __ _ __  __Q_QA___-EEAá<x<xx==5 
onto a óqxqq- (_ _ _— ——_———————— 
€ ___z_—_—___z---_ < «<__—_  —_AAA=>> E (j]h__EOE EE EAá- A A A A 
(A e A _ — ___ __ 
¡Egg zm _xí_ _—_ ___-_ —_—_— ———_ — —  _— ————— 
a _-z__—__ ___——_—_—__ _—_——__—_ —_—__—_ _——_—_— _ _—_ —__————- ¿€ _z A A A 
O ___ ———_ _ __=-z-M« _ __—_ __ —_—_ _ — — _ _  _——-- _  ___  _---==44 
nee ——Á 
x€TOTOnB BR _z-__ o_o _—_—— _—_— _——— — — _ —— —_——_—  _ _ __ _>_<-xH3 HH 
y ]__—a_awm=z=m=wxwmxRm_  _ __ mm —___—____—___—___—_ ___—__ __ _ Q_»>>oa_ 
A A q 
e A _ _____ ___-_-__--_——_ KA 
q _OÓÉ]AÓ_ooeoeoo Qr. .IRooF=oPFe _ -_- xh____ a —_—_  _—_ —___ _ __—_—__-_»->___ A A A 
€ _-__-_-__ A  __ _—_  ___________ __ A A A AAA 
e e e xq í____— _— A 
a  __ _____ _ — ___—  _—__— _—___—_ _ —___—__ _A>E_ÁA __--=>z=>——2— 


e Defina las matrices: Las variables definidas en PASCAL como ARRAY 
debe definirlas en BASIC con una sentencia DIM: 


e Convierta las funciones de entrada/salida: 


e WRITE o WRITELN deberá transformarla por PRINT. 
+ READ o READLN deberá transformarla por INPUT. 


* Sentencias de bucles (FOR) son muy fáciles de trasladar: 
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REPEAT 10 sentencia 


sentencia; 20 sentencia 
sentencia; 30 sentencia 
sentencia; 40 IF C$<>"Q"” THEN GOTO 10 


UNTIL Componente='0Q' 


WHILE Control <> 'N' do 10 IF C$ = "N" THEN GOTO 60 
begin 20 sentencia 
sentencia; 30 sentencia 
sentencia; 40 sentencia 
sentencia; 50 GOTO 10 
end; 60 sentencia 


e Procedimientos. Esta es la parte más delicada, ya que el BASIC no 
dispone de llamada a subrutinas con paso de parámetros. En cualquier 
caso, hay que dar a las variables que se van a utilizar en el procedimiento 
los valores para los que queremos que dicho procedimiento se ejecute. 


Veámoslo de forma más clara mediante un sencillo ejemplo: 


PROGRAM Demostracion; 


VAR 
m,n,o: integer; 


PROCEDURE Ejemplo (VAR a,b,c: INTEGER); 
VAR 1: integer; 
BEGIN 

1:=b*c+3; 

IF 1<0 THEN 


BEGIN 
READLN(m,n,o); 
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A A 
lempioj! MM, 07) 


URTTELN(mM 1O7J 
¿AA Er 


| 


Es muy importante resaltar que las variables definidas dentro de un pro- 
cedimiento en PASCAL (la «1» del programa anterior, por ejemplo), son va- 
riables locales. Esto quiere decir que no son visibles desde el resto del pro- 
grama. Por tanto: 


e Su valor sólo es accesible dentro del procedimiento. 

e Si existiera otra variable con igual nombre en el resto del programa, 
su valor no se vería alterado al asignarle un valor a la variable local en 
cuestión. 


En BASIC no existen las variables globales. Cuide, por tanto, su uso den- 
tro de los procedimientos. Si en el ejemplo anterior «l» fuera una variable 
del programa principal, su valor se perdería debido a la sentencia de asig 
nación de la línea 1000. 

A continuación se ofrece el programa de resolución de sistemas de 
ecuaciones convertido al BASIC como ejemplo: 


LET OM = 20 
REM Matriz del sistema 
DIM A(0M,0M) 


REM Vector de términos independientes. 
DIM B(0M) 
PRINT "RESOLUCION DE SISTEMAS DE ECUACIONES” 
PRINT Un ----- = 
PRINT 
PRINT "Orden del sistema...... == 
INPUT N 
PRINT 
PRINT "Introduzca los coeficientes del sistema:” 


) FOR I=1 TO N 
FOR J=1 TO N 
PRINT ACI AI As > 
INPUT A(1,J) 
NEXT J 
NEXT I 
PRINT; 
PRINT "Introduzca los términos independientes:” 
FOR I=1 TO N 
ERRE == 
INPUT B(1) 
NEXT 1 
GOSUB 1000 
PRINT'; 
PRINT("'Descomposición del sistema:”); 


PRINT'; 

PRINT("Matriz L:”) 

FOR-I=1-TO-=N 

FOR J=1 TO I 

PRINT AZ = A 
NEXT J 

NEXT 1 

PRINT 

PRINT "Matriz U:" 


FOR I=1 TO N 

FOR J=I+1 TO N 

PRINT SEAEESAES = "sA(1,3) 
NEXT J 

NEXT I 

PRINT 

PRINT "Solución del sistema:"” 


FOR I=1 TO N 


Iposicion A = Ll 
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La simulación por ordenador en 
general, y la de circuitos en particular, 
es uno de los temas que mayor interés 
han despertado en los técnicos desde el 
comienzo de la informática hasta 
nuestros días. En el presente libro se 
estudian los métodos matemáticos que 
permiten representar un circuito 
mediante fórmulas y ecuaciones, y las 
herramientas necesarias para 
aprovechar al máximo la potencia del 
ordenador con vistas a automatizar el 
proceso de diseño. 


